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Abstract. We consider the problem of quantization of classical soliton integrable 
systems, such as the KdV hierarchy, in the framework of a general formalism of 
Gaudin models associated to affine Kac-Moody algebras. Our experience with the 
Gaudin models associated to finite-dimensional simple Lie algebras suggests that the 
common eigenvalues of the mutually commuting quantum Hamiltonians in a model 
associated to an affine algebra g should be encoded by afhne opers associated to the 
Langlands dual affine algebra ^g. This leads us to some concrete predictions for the 
spectra of the quantum Ifamiltonians of the soliton systems. In particular, for the 
KdV system the corresponding affine opers may be expressed as Schrodinger opera- 
tors with spectral parameter, and our predictions in this case match those recently 
made by Bazhanov, Lukyanov and Zamolodchikov. This suggests that this and other 
recently found examples of the correspondence between quantum integrals of motion 
and differential operators may be viewed as special cases of the Langlands duality. 



1. Introduction 

Soliton equations, such as the celebrated KdV hierarchy, are infinite-dimensional 
classical integrable systems possessing infinite collections of Poisson commuting Hamil- 
tonians. It is natural to try to quantize these systems. This means, in particular, 
constructing mutually commuting quantum Hamiltonians in the appropriate associa- 
tive algebra (such as the completed enveloping algebra of the Virasoro algebra in the 
case of the KdV hierarchy), whose symbols are the classical Hamiltonians, and finding 
the spectra of these quantum Hamiltonians on representations of this algebra. This 
problem has acquired additional importance after the work of A. Zamolodchikov [Z], 
in which he showed that these quantum Hamiltonians may be identified with integrals 
of motion of certain deformations of conformal field theories. 

The existence of local commuting quantum Hamiltonians has been established in 
|FF41 IFF5j in the case of the generalized KdV hierarchies and in [FF6j in the case 
of the generalized AKNS hierarchies. The corresponding non-local Hamiltonians have 
been constructed in [ BLZlt IBHKj for the KdV hierarchies (we will discuss below an 
extension of this construction to the AKNS hierarchies). A study of the problem of 
diagonalization of the quantum KdV Hamiltonians was initiated in the important series 
of works |BLZ1| - |BLZ3| by V. Bazhanov, S. Lukyanov and A. Zamolodchikov. 
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These works culminated in a remarkable conjecture made in |BLZ5| (building on the 
earlier works [DTll lBLZ4j ). which related the spectra of the quantum KdV Hamiltoni- 
ans to certain one-dimensional Schrddinger operators. This connection appears to be 
rather mysterious. In this paper we suggest that it may be viewed as a special case of 
a broader picture which fits in the general framework of integrable systems of Gaudin 
type and the Langlands correspondence. 

The first step in our approach is the realization that the quantum soliton systems 
may be viewed as special cases of Gaudin models associated to affine Kac-Moody al- 
gebras. At the classical level, this is quite clear. Indeed, it has been known since 
[RSH IDS| that the algebra of classical Poisson commuting Hamiltonians of the soli- 
ton hierarchies may be obtained using the "shift of argument" method, applied to an 
affine Kac-Moody algebra. Its finite-dimensional counterpart is the Poisson commuting 
"shift of argument" subalgebra constructed in [MF] . As explained in [FFTj , the corre- 
sponding Hamiltonian system is the simplest example of a classical Gaudin model with 
irregular singularities (which is in turn a special case of the Hitchin systems jiH l IBea| ) . 
Recently, a quantization of this algebra, as well as its generalizations corresponding to 
other Gaudin models, has been constructed in [R[ [FFTj. We expect that the soliton 
integrable systems may be quantized along the same lines, that is, in the framework of 
the quantization of Gaudin models. 

The second step in our approach is the link between the spectra of the quantum 
Gaudin Hamiltonians and differential operators known as opers, discovered in our ear- 
lier works [FFRl \F3\ IFFTj . We have shown there that the spectra of the Gaudin 
Hamiltonians associated to a simple finite-dimensional Lie algebra q are encoded by 
^g-opers (where is the Langlands dual Lie algebra of g) on with finitely many 
singular points. 

The key point of our analysis |FFRl IF3t IFFTj of the Gaudin models associated to a 
simple Lie algebra g is the identification of the commutative algebra of quantum Gaudin 
Hamiltonians with a quotient of the center Z{q) of the completed enveloping algebra 
of the affine Kac-Moody algebra g at the critical level. Now we consider the Gaudin 
systems associated to an affine algebra g. Therefore it is natural to assume that the 
corresponding Hamiltonians come from the center Z(q) of the completed enveloping 

algebra of the double affine algebra g. At present, we do not have a precise definition of 
this enveloping algebra or its center. However, we postulate that it exhibits the same 
salient features as its affine counterpart Z(g). 

According to a theorem of |FF3tlF4] . Z(q) is isomorphic to the algebra of functions on 
the space OpLg(D^) of ^g-opers on the (formal) punctured disc . Hence we expect 

that Z{g) is related to the space of affine opers OpLg(D^) associated to ^g, defined in 
|F2| (see also [BFj). Here ■'"g is the Langlands dual Lie algebra to the affine algebra 
g, that is, its Cartan matrix is the transpose of that of g (so the dual of an untwisted 
affine algebra might be a twisted affine algebra). Testing this conjectural relation was 
in fact one of the main motivations for the present work. Since the structure of the 
center Z(q) seems out of reach at the moment, we turn the tables and try to use the 
available information about the affine Gaudin models and the quantum soliton systems 
to gain insights into the structure of Z(q). 
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The upshot of all this is that the spectra of the quantum g-KdV Hamiltonians should 
be represented by the affine ^g-opers on with prescribed singularities. In the case 
when Q = 5l2 these opers may be represented by Schrodinger operators with a spectral 
parameter. Thus, we suggest that the connection between the spectra of quantum KdV 
Hamiltonians and Schrodinger operators conjectured in |BLZ5] comes about as follows: 

• The quantum KdV system is interpreted as a generalized Gaudin model asso- 
ciated to the affine Kac-Moody algebra 5(2 (or g for the more general quantum 
g-KdV systems). 

• The spectra of the Gaudin Hamiltonians associated to an affine Kac-Moody 
algebra g are encoded by the affine ^g-opers. 

• Since the quantum KdV Hamiltonians may be interpreted as the Gaudin 
Hamiltonians associated to SI2, their spectra should be encoded by s[2-opers, 
which are nothing but Schrodinger operators with spectral parameter: 

— v(t) — A. More generally, we expect that the spectra of the quantum 
g-KdV Hamiltonians are encoded by affine ^g-opers. These opers should sat- 
isfy an important no monodromy property. 

Thus, schematically, the correspondence we propose looks as follows: 



quantum KdV 
Hamiltonians 



affine Gaudin 
models 



opers for the Langlands 
dual affine Lie algebra 



The appearance of the Langlands dual Lie algebra suggests that this connection 
(and perhaps the more general ODE/IM correspondence of Dorey, e.a., [DDTHlDMSTl 
IDDT2j ) should be viewed as a special case of Langlands duality. 

In order to test this proposal, we compare our 5l2-opers with the Schrodinger op- 
erators obtained in [DTH IBLZ4[ IBLZ5] . We find that they match perfectly, after a 
simple change of variable. Our formulas easily generalize to an arbitrary affine Lie 
algebra g, and again we find that they are consistent with the known formulas for 
the differential operators associated to the spectra of the quantum Hamiltonians in 
[DDTli "BHIT, ' DMSTIIDDT2]. 

In [DTI. .BLZil lDT2l IDDTll iBHKl IBLZ5j the connection between the spectra of 
the quantum KdV Hamiltonians and differential operators is explained by the fact 
that certain generating functions of the former (the so-called Q-operators) satisfy the 
same functional relations as connection coefficients (or Stokes multipliers) of the latter. 
These functional relations are the so-called T — Q relations generalizing Baxter's famous 
relation. In |DMSTllDDT2j (and earlier works referenced there) the functional relations 
of this type were taken as the basis of an "ODE/IM correspondence" for more general 
solvable models. 

In contrast, in this paper we do not attempt to explain the connection between the 
spectra of quantum Hamiltonians and affine opers. Rather, our goal is to emphasize 
that the connection between them fits in the general scheme of quantum Gaudin models 
and is in fact a manifestation of the Langlands duality. This is based on our analysis 
[FFR1IF31IFFT] in the finite-dimensional case and on the fact that in the affine case our 
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predictions match those obtained in [DTT| IBLZ41 IDDTH \BBK\ IBLZ5j . This already 
has a number of consequences: 

(i) the appearance of the Langlands dual affine Lie algebra (something that would 
be difficult to see based on the most studied case of s[„, which is self-dual); 

(ii) a uniform way of expressing the differential operators encoding the spectra of 
the quantum Hamiltonians as opers associated to arbitrary affine Lie algebras 
(again, this is not obvious from the previously known examples); 

(iii) a Bethe Ansatz method for constructing eigenvectors of the quantum Hamil- 
tonians. 

In particular, we conjecture explicit Bethe Ansatz formulas (generalizing the Bethe 
Ansatz in the finite-dimensional Gaudin models) for eigenvectors of the quantum KdV 
Hamiltonians in the irreducible representations of the Virasoro algebra from the unitary 
minimal models. 

We expect that the appearance of the T - Q relations of [DTT| IBLZ41 [DT2l IDDTll 
IBHKl IBLZ5| may also be explained in the framework of the formalism developed in 
this paper. More specifically, we believe that it comes from the connection between 
quantum integrals of motion and the center of a quantum affine algebra at the critical 
level. We plan to address these issues in a future work. 

The paper is organized as follows. In Section [2] we recall the results of [R l IFFT] on 
the quantization of the shift of argument Gaudin algebras in the finite-dimensional case 
and the identification of their spectra with ^g-opers. We then consider in Section [3] 
the affine analogue of this construction which corresponds to the quantization of the 
Hamiltonians of the soliton hierarchies such as the g-AKNS. We discuss in detail the 
classical Hamiltonians obtained by expansion of the monodromy matrix with respect to 
the spectral parameter (following [RSlJ) and the quantum analogues of both local and 
non-local Hamiltonians. We then construct eigenvectors of these Hamiltonians using 
an affine analogue of the Bethe Ansatz method. In order to describe the spectra of the 
quantum Hamiltonians, we introduce in Section H] the affine opers, following [F2j. We 
conjecture that the affine opers that encode the spectra of the quantum Hamiltonians 
on representations of g are ^g-opers on satisfying certain analytic properties (in 
particular, to integrable representations of g correspond ■^g-opers without monodromy). 

Next, in Section [H we apply the same analysis to the generalized KdV hierarchies. 
We state our conjectures describing the spectra of the quantum KdV Hamiltonians in 
terms of the affine opers corresponding to the Langlands dual affine Lie algebra ^g. We 
write down explicit formulas for these opers in the case when g = s[2 and show that after 
a simple change of variables they give rise to the Schrodinger operators constructed in 
[BLZ5j . We then discuss the generalization of these formulas to the case of s[„ and to 
general affine Lie algebras. 

Motivated by this study of the "shift of argument" soliton systems, corresponding to 
the simplest affine Gaudin models with irregular singularities, we introduce in Section [6] 
the affine analogues of the Gaudin models with regular singularities. We show that in 
this case we also have interesting classical integrable systems, which may be quantized 
by analogy with the finite-dimensional case. Again, we conjecture that the spectra of 
the corresponding quantum Hamiltonians are described by the affine ^g-opers on 
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with regular singularities at finitely many points. We also construct the eigenvectors of 
these Hamiltonians using Bethe Ansatz. We show that a particular case of this model is 
closely related to the GKO coset construction |GKO) . We conjecture that in this case 
our Bethe vectors are eigenvectors of all quantum KdV Hamiltonians on irreducible 
modules over the Virasoro algebra from the unitary minimal models. 
We summarize our results and conjectures in Section [71 

Acknowledgments. This project was conceived and to a large extent carried out 
during our visits to Kyoto University in 2004-2007. We thank our hosts M. Kashiwara, 
T. Miwa and H. Nakajima, as well as K. Hasegawa, for their hospitality. 
We also thank V. Bazhanov for useful discussions. 

2. The shift of argument Gaudin models in the finite-dimensional gase 

We start by summarizing what is known about the simplest Gaudin models with 
irregular singularities, following [R[ IFFT| . The corresponding classical Hamiltonian 
system may be viewed as an example of a Hitchin system. The Poisson commutative 
algebra of classical Hamiltonians is the so-called "shift of argument" subalgebra Ay. of 
the symmetric algebra S{q) of a simple Lie algebra g, introduced by A.S. Mishchenko 
and A.T. Fomenko |MFj (see also [Ml]). The corresponding algebra of quantum Hamil- 
tonians is a commutative subalgebra of the universal enveloping algebra U{q). It has 
been constructed in [R[ IFFTj using the center of the completed universal enveloping 
algebra of the affine Kac-Moody algebra q at the critical level. This allowed us to iden- 
tify in [FFTj the spectrum of this algebra with the space of ^g-opers on with regular 
singularity at € and irregular singularity at oo G P^ (here is the Langlands 
dual Lie algebra to g). 

We also consider here analogous commutative subalgebras in the Hamiltonian reduc- 
tions of C/(g) known as finite W-algebras. 

In the following sections we will use this construction as a prototype for the de- 
scription of the quantum "shift of argument" subalgebra when the finite-dimensional 
Lie algebra g is replaced by its affinization g. It is well-known (see jRSH IDSj ) that 
the corresponding classical systems give rise to hierarchies of soliton equations such as 
g-AKNS and g-KdV. In other words, the classical shift of argument subalgebra is gen- 
erated by the classical Hamiltonians of these hierarchies. Therefore the corresponding 
quantum Hamiltonians generate the quantum shift of argument subalgebra, and this 
suggests that their spectra are encoded by affine ^g-opers. We will find that in the case 
of the quantum KdV Hamiltonians this matches the predictions of |BLZ5| . 

2.1. Classical system. The Gaudin models with irregular singularities for finite- 
dimensional simple Lie algebras have been introduced and studied in |FFT] . Here 
we will focus on the particular case related to the so-called shift of argument subal- 
gebra A-)^ C 5(g) and its quantization A^ C U{q), as it is most closely related to the 
structures emerging in soliton equations (such as the KdV and n-wave hierarchies) and 
their quantization, once we switch from g to the affine algebra g. But other Gaudin 
models with irregular singularities may also have interesting analogues in the affine 
case. 
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Our starting point is a classical Hitchin system (see [H| IBeaj ) in which the phase 
space is the space of the following g*-valued one- forms on P^: 

(2.1) i]=(^^+x^dz, A eg*, 

where x ^ 0* is fixed. This space is therefore isomorphic to g*. A one-form like this 
has a pole of order 2 at the point oo with the leading coefficient —x, and it also has a 
pole of oder 1 at the point 0. 

We define the classical shift of argument algebra as the subalgebra of S{q) = 
Fung* generated by the coefficients of the polynomials in obtained by evaluating 
the invariant polynomials P € (Fun Q*)^,i = 1, . . . ,i, on rj. In other words, is the 
subalgebra of Fung* generated by the iterated directional derivatives D^P of invariant 
polynomials P E (Fung*)*^ in the direction x, where 

d 



(2.2) D,P(,) 



P{y + ux)- 

u=0 



Equivalently, A^ is the subalgebra of Sq generated by the shifted polynomials 

(2.3) Pu^{y) = P{y + ux) 

where P varies in (Fung*)*-' and ti G C, hence its name. 

The algebra A^ is Poisson commutative with respect to the Kirillov-Kostant Poisson 
structure on g*. This follows from the general results on the Hitchin systems with 
singularities. 

The algebra A^ was first introduced in [MFj . It has been shown in [MF^ IKo^ IFFT) 
that for regular x the algebra A^ is a maximal Poisson commutative subalgebra of S{q) 
which is a free polynomial algebra in dim b generators 

(2.4) D'^^Pi, i = l,...,e; n^ = 0,...,di, 

where Pi is a generator of S(q)^ of degree di + 1 and is the derivative in the direction 
of X given by formula (j2.2p Pl 

Note that for any x the algebra A^ belongs to the centralizer S{q)^^ of in S'(g), 
where g^^. C g is the stabilizer of x- Let ip G g*^he a character of g^^., or, in other words, 
a one-point coadjoint orbit in g* . Set 

5^(g) = (5(g)/(g^-V(gx)))^^ 
where (g^^ — V'(Sx)) ^^e ideal in S{q) generated by the elements A — ip{A),A G g^. 
The algebra ^^f (g) is a Poisson algebra obtained by Hamiltonian reduction of g* with 
respect to the adjoint action of Gy. and the one-point orbit ip G g*. Let A^ be the 

projection of A^^ onto S^{q). This is a Poisson commutative subalgebra of 5'^(g). 

For regular x the Lie algebra g^^ is commutative, and its symmetric algebra S{q^) 
is contained in Ay.. In fact, a basis of g^^- may be obtained by taking the maximal 
non-constant directional derivatives W^'Piji = !,...,£, among the generators (12. 4p . 



"'^We recall that a £ g is called regular if its centralizer in g has the smallest possible dimension; 
namely, £, the rank of g. We identify g* with q by using a non-degenerate invariant inner product. 
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Therefore for any i/j G g* the algebra Ay, is a free polynomial algebra generated by the 



(non-zero) images of 



D^'Pi, i = l,...,e; ni = 0,...,d, 



in St (9). 

For nilpotent elements x the algebras 5*^(0) are known as the (classical) finite W- 
algebras, and the above construction produces Poisson commutative subalgebras in 
these W-algebras. 

For example, suppose that x corresponds to the maximal root generator G n_ C g 
under an identification g ~ g* given by a non-degenerate invariant inner product on 
g. Then g-^ = n_. Let ■0 be a principal character on n_ taking non-zero values on the 
simple root generators /j. We have ^^(g) = 'S'(g)^ in this case, and therefore it is clear 
that = ^(g)^ as well. However, the affine analogue of this construction is non-trivial: 
the analogue of ^^(g) is the classical W-algebra obtained by the Drinfeld-Sokolov 

reduction of g^, and the analogue of is its Poisson commutative subalgebra generated 
by the (local and non-local) Hamiltonians of the g-KdV hierarchy (see Section E]). 

2.2. Quantization. Let us now discuss the quantization of the algebra A^ in the finite- 
dimensional case. Here by quantization we understand a commutative subalgebra Ay. of 
U (g) such that the associated graded gr ^1;^ of A^ (with respect to the standard filtration 
on C/(g)) is a Poisson commutative subalgebra of grC7(g) = S{q) which contains Ay.. 
Such a quantization has been constructed in [R[ IFFT] in a uniform way for an arbitrary 
simple Lie algebra g (following the procedure that had been used for constructing the 
algebra of quantum Hamiltonians in the Gaudin models with regular singularities in 
[FFR1IF3] . see Section O below) 

Let g be the affine Kac-Moody algebra associated to g (see Section [3T|) . and 

(2.5) Vcnt = Ind;[[,jjg^^,C.,v 

the vacuum g-module of critical level k = —hy, where is the dual Coxeter number 
of g (in this formula g[[t]] acts by on the one-dimensional module C_/iV, and 1 acts 
by multiplication by — /i^). The quantum shift of argument algebra Ay. was obtained in 
[R1 IFFT] as a quotient of the algebra End^ Vcrit- According to a theorem of |FF31 |F4] . 
the latter algebra is isomorphic to the algebra of functions on the space OpLg(-D) of 
^g-opers on the (formal) disc. Here ^g is the Langlands dual Lie algebra to g, whose 
Cartan matrix is the transpose of the Cartan matrix of g. The notion of oper was 
introduced by Beilinson and Drinfeld [BDll [BD2| following Drinfeld and Sokolov [DSj . 
We refer the reader to the definition given in ]F2\ IF3j . which is well-suited for our 
present goals. (We will recall this definition below in Section SI) 
Thus, we have an isomorphism 

(2.6) End^Vcrit ^ FunOpL„(L'). 



o — 

For classical Lie algebras the quantization of had been constructed previously in a different way 
in [NOimiCT]. 
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We note that End^ Vcrit is the quotient of the center Z(q) of the completed enveloping 
algebra of g at the critical level, which is isomorphic to the algebra of functions on the 
space OpLg{D^) of opers on the (formal) punctured disc (see, e.g., |F6] . Sect. 4.3). 

Using the isomorphism (j2.6p and general results on coinvariants from [FBJ , we have 
obtained the following description of Ay. in |FFTj : 

Theorem 1. For regular x G 0* the algebra is canonically isomorphic to the algebra 
of functions on the space OpLg(P"'^)^(^) of ^g-opers on with regular singularity at the 
point and with irregular singularity of order 2 at oo, with the 2-residue ir{x)-' 

(2.7) A^c^OpL^{FX(^yy 

For example, in the case when g = 5(2 and x £ f)*) the space OpLg(P^)^(^) consists 
of second order differential operators of the form 

(2.8) a2_^2^!^^!^^ v_,,v_,eC. 

Note that they have a pole of order 4 at z = cxo with the leading term — x^- For g = st^ 
the role of operators ()2.8p is played by the following nth order differential operators: 

where Vj{z) has the form 

i+i 

Vj{z) = Pj{x) + ^Vj^kZ~^, 

k=l 

and Pj is a particular symmetric polynomial of degree j + 1- 

For an explicit description of OpLg(P^)^(^) for general simple Lie algebras, see |FFT] . 

It has been proved in [R l IFFT ] that gr^lj^. contains A^, with the equality for regular 
X- For general x ^ 9* this equality is a conjecture (see Conjecture 1 of |FFTj). 

Let us summarize: the quantization A^ of the classical Poisson commutative algebra 
A^ is a commutative subalgebra in a quantum (non-commutative) algebra U{q), but 
it turns out to be isomorphic to another classical commutative algebra, namely the 
algebra of functions on the space OpLg{F^)T^(^yy A conceptual understanding of this 
isomorphism is obtained by realizing Ay. as a quotient of the center of an enveloping 
algebra of the affine Kac-Moody algebra g and using the identification of the latter with 
the algebra of functions on opers. The non-triviality of this isomorphism is underscored 
by the fact that the Lie algebra g gets replaced by its Langlands dual Lie algebra ^g 
under this isomorphism. In fact, this isomorphism may be viewed as an example of the 
geometric Langlands correspondence, as explained in [Fit IF5t IFFT] . 

The Poisson commutative subalgebras A^ C ^^(g) may also be quantized. Namely, 

by construction, Ay^ is contained in U{q)'^^. For a character -0 : g^ ^ C, let Ay be the 
projection of Ay onto 

^x^(0) = (c^(0)/(0x-V'(0x)))'^ 

where (0^ — ip{5y)) is the left ideal in U{g) generated hy A — 'ip{A),A € Qy. Thus, 
Uy{g) is the result of quantum Hamiltonian (or BRST) reduction of U{g) with respect 
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to Qy^ and its character ip. In particular, if x is nilpotent, U^{q) is a quantum finite 
W-algebra, and so Ay^ gives rise to its commutative subalgebra. For example, if x 
corresponds to the maximal root generator /e € n_ C and ■0 is a principal character 
on n_ (as in Section [2?T]) . then = U^ig) = Z{q), the center of U{g). 

2.3. The spectra of the quantum Hamiltonians. Suppose that X ^ 0* — is 
regular semi-simple. Without loss of generality we may assume that x is in a fixed 
Cartan subalgebra f) C 0. The algebra Ay. then contains f). It is natural to try to 
describe the spectrum of A^ on a given g-module M. Note that since i) C A^, we obtain 
that the action of A^. will preserve the weight decomposition of M. The isomorphism 
(j2.7p implies that the common eigenvalues of A^. on M are encoded by ^g-opers on 
with regular singularity at the point and with singularity of order 2 at oo, with the 
2-residue tt{x) (which is, by definition, an element of i)* /W, see, e.g., |FFT) ). In |FFT] . 
Theorem 5.10, we have obtained the following more precise description of the common 
spectra of A^. on the Verma modules and irreducible finite-dimensional representations 
of 0. 

Theorem 2. 

(1) If M = My, the Verma module with highest weight v ^ \)* , then the common 
eigenvalues of A^ on M are encoded by ^g-opers on with singularity of order 2 at 
oo, with the 2-residue tt{x) and with regular singularity at the point and the residue 
(which is, by definition, an element of \)* /W , see, e.g., \F YT\ ) equal to w{—i' — p), 
where w is the projection f)* ^ \)* /W . 

(2) If M = Vy is an irreducible finite- dimensional representation of q, then the 
common eigenvalues of A^ on M are encoded by the ^Q-opers in part (1) satisfying the 
additional condition that they have trivial monodromy. 

According to Conjecture 2 of [FFTj . all opers described in part (2) of this theorem 
are realized in the spectrum of A^ on Fj^H 

This is analogous to the description of the spectra of the Hamiltonians in the Gaudin 
models with regular singularities from |FFR1 lF3j which we will recall in Section 16.11 

In [FFTj we have constructed eigenvectors of Ay. by Bethe Ansatz. Suppose that 
M = My. Then the Bethe eigenvectors have the form 

(2.9) </.(<,..., «;j;r) = 

/jo-(l)/*CT(2) • • • f'i-cr(m) 



E 



where Vy is the highest weight vector of My, and the sum is over all permutations on 
m letters. This vector has the weight 



This conjecture has been proved by L. Rybnikov and the authors in Opers with irregular singularity 
and spectra of the shift of argument subalgebra, Preprint arXiv:0712.1183 
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The following result is proved in [FFT] . Sect. 6.6 (for the definition of Miura trans- 
formation, see |FFT] and Section [4] below). 

Theorem 3. The Bethe vector 4>{w^i , ■ ■ ■ , w^^) is an eigenvector of if the complex 
numbers Wj satisfy the Bethe Ansatz equationi^ 

(2.10) + ^^)=0, J = l,...,m. 



Wj — Ws 



The corresponding eigenvalues of elements of Ay, are encoded by the ^g-opers in 
OpLg(P"'^)^(^) obtained by applying the Miura transformation to the connection 



m 



(2.11) v = a,-x-- + 5^ 



Z ' Z — Wn 



on the ^H-bundle on 



ii 



Conjecturally, these Bethe eigenvectors form an eigenbasis of Miy for generic x 
z^, and a basis of for a dominant integral i/ and generic x (iii this case the above 
weights should be in the set of weights of Vi,). 

Finally, we note that some of the Hamiltonians in A^ may be constructed explicitly. 
Assume again that x ^ f) C g is regular semi-simple. In this case the span of all 
elements of Ay- of order less than or equal to 1 (with respect to the standard filtration 
on U{q)) is the direct sum of the Cartan subalgebra = Qx 9 U{q) and the scalars 
C C U{q). Next, we have the following quadratic elements: 

(2.12) T,{x) = Yl 7^ + ' 7 G r , 

where the elements Ca, fa generate the 5(2 subalgebra Qa of q corresponding to a G A+, 
and they are normalized in such a way that KQ^Ca, fa) = 1; where kq is a non-degenerate 
invariant inner product on g. These are the DMT Hamiltonians introduced by C. De 
Concini (unpublished) and by J. Millson and V. Toledano Laredo |MTH ITL] (see also 
[FMTV ])PI Explicit formulas for higher order Hamiltonians are unknown in general. 
However, for 0=5!^ explicit formulas have been found in [T j ICT| . 

3. The shift of argument affine Gaudin models 

We now try to generalize the above construction to the situation where a finite- 
dimensional simple Lie algebra g is replaced by the corresponding (extended) affine 
Kac-Moody algebra g. 

In what follows, we will focus on the case when g is an untwisted affine algebra. 
However, all of our definitions, results and conjectures presented below generalize in a 
straightforward way to the case of twisted affine algebras. 



^Note that our x here corresponds to — x in [FFT) . 

^We note that quasi-classical versions of these Hamiltonians have appeared in [Boa] in the study of 
isomonodromic equations with irregular singularities. 
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3.1. Affine Kac— Moody algebras. Let 

(3.1) ^ CI ^ ^ g{{t)) ^ 

be the Kac-Moody central extension of Q{{t)) with the commutation relations 

[A f{t),B ® g{t)] = [A, B]®fg- ko{A, B) Rest=o fdg ■ 1, 
=0, 

and kq is the invariant inner product on g normalized in such a way that the maximal 
root has squared length 2. (We will sometimes write {A,B) instead of kq{A,B).) Let 
g be the corresponding extended affine Kac-Moody algebra 

= Cd x 0, 

where d acts trivially on 1 and as the outer derivation tdt on g((t)): 

[d,A0fit)]=A(^tdtf{t). 
Let us fix a Cartan decomposition 

s = n_ e [) en+. 
Then the Lie algebra q has the following Cartan decomposition: 

= n_ e^en+, 

where 

(n+ 1) e {g(^tC[[t]]), 
(n_ (g) 1) e {g0t-^C[t-^]), 

(i) o 1) ecl ecd 

The Lie algebra g is a symmetrizable Kac-Moody algebra (see [K]). We denote by / 
the set of vertices of the Dynkin diagram of g. Thus, I = I U {0}, where I = {1, . . . , i} 
is the set of vertices of the Dynkin diagram of g. We have the generators of n+, fi 
of n_, /ij = Oi and 1 of f), where in all cases i £ I. We also denote by a, G f)*, i G /, 
the simple roots of g. 

A generalization of the shift of argument method to affine Kac-Moody algebras 
is well-known to give rise, classically, to Poisson commuting Hamiltonians of a large 
class of soliton equations. This point of view on soliton equations was explained by 
A. Reyman and M. Semenov-Tian-Shansky ^RSlj (see also related works {M2t IRSF] ) 
and has been widely used ever since. In particular, this construction was subsequently 
generalized by V. Drinfeld and V. Sokolov |DS] to include Hamiltonian reductions 
(this yields the generalized KdV hierarchies, as we explain in Section [5] below). In 
this section and the next we will discuss the affine Gaudin models, both classical and 
quantum, corresponding to the regular semi-simple shift. They realize the g-AKNS 
soliton models (or n-wave systems) . Then in Section [5] we will consider the g-KdV 
models corresponding to the Drinfeld-Sokolov reduction of the shift by the nilpotent 
maximal root element. 



n+ = 
n_ = 
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3.2. The classical model. Consider the case when the shift is a regular semi-simple 
element 

X G f) C C g. 

The phase space is then an affine analogue of the space of one- forms (j2.ip , so it consists 
of the one-forms 

(3.2) r^=(^^+^^dz, A eg*- 

Instead of the full dual space q* will consider its codimension 2 affine subspace which 
consists of the functionals taking the value k ^ on the central element 1 and on the 
element d. In other words, this is the affine hyperplane C g* which consists of the 
linear functionals on g taking value A; on 1. It is endowed with the Kirillov-Kostant 
Poisson structure, and it is well-known that it may be identified with the space of 
fe-connections 

(3.3) kdt + A{t), A(t)Gg*((t))^g((t)) 

(where the second isomorphism is obtained using the non-degenerate invariant inner 
product kq on g) on the punctured disc, so that the adjoint action of G{{t)) on 
becomes the action of G{{t)) on fc-connections by gauge transformations (see [RSHlFB] ). 

By rescaling A{t) i— > ^A{t), we may identify g^ and gj, so without loss of generality 
we may, and will, set = 1. 

Thus, we replace rj in (|3.2p by 



and further by the equivalent expression 

(3.4) L = dt + A{t) + xz. 

We now recognize it as the L-operator of the so-called the g-AKNS hierarchy, also 
known as the generalized n-wave hierarchy (see, e.g., [RSH [D5] ). 

The Poisson commuting Hamiltonians of the corresponding integrable system are 
obtained by expanding in z (or z~^) the invariants of the monodromy of the connection. 
Denote by M{z) € G (where G is the connected simply-connected Lie group with the 
Lie algebra g) the formal monodromy matrix of the connection operator L given by 
formula ()3.4p . More precisely, M{z) is defined as follows: its action on any finite- 
dimensional representation y of G is obtained by computing the monodromy of the 
matrix differential operator dt + A{t) + xz acting on y|j 

Now, for any function 99 on G which is invariant under the adjoint action of G, the 
Hamiltonian 

H^{z) = ip{M{z)) 



*^Strictly speaking, this monodromy is defined for connections that converge on a small punctured 
disc around 0; however, the coefficients in the expansions that we will consider are also well-defined 
when A{t) is a formal Laurent power series. 
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is a gauge invariant function of L given by formula (13. 4p . depending on the spectral 
parameter z. One shows in the same way as for the Hitchin systems in the finite- 
dimensional case that 

(3.5) {H^{z),H^{y)} = {), V(/p, V e Fun(G)G, z, y G C\ 

Therefore we may obtain infinite collections of commuting Hamiltonians of the affine 
Hitchin model by expanding the Hamiltonians H^{z) around the points and oo € P^. 

The asymptotic expansion of H^{z) at z = oo gives rise to the local Hamiltonians, 
in the sense that they may be expressed as integrals of differential polynomials in the 
matrix coefficients of A{t) (in the adjoint representation). There are standard methods 
for constructing these Hamiltonians explicitly, see [RSH IDS! IFF6j . These Hamiltonians 
generate the generalized g-AKNS (or n-wave) hierarchy. Since the algebra of invariant 
functions on G is a polynomial algebra with i = rank(g) generators, we find that there 
are t local Hamiltonians of each positive spin. 

On the other hand, the expansion of H^{z) around z = gives rise to the non-local 
Hamiltonians. Let us consider the simplest examples of such Hamiltonians obtained by 
taking the z-linear terms of H^{z). 

The computation of these linear terms is based on the following observation which is 
well-known in the study of the inverse scattering problem (for a proof, see, e.g., |RS2] . 
Lemma 4.2.4). 

Lemma 1. Suppose that we have a connection of the form 

(3.6) 5^ + A(0)(t) + eA«(t) + ..., 

Then the e-linear term in the expansion of ip{M {e)) , where M(e) £ G is the monodromy 
of this connection and (p is an invariant function on G, is equal to 

(3.7) l{AW{t),^^{t))dt 

(the integral is taken over a contour around 0), where ^'<^(t) is a solution of the equation 

(3.8) 9t*^(t) + [AW(t)),^^(t)]=0, 

which is single-valued (i.e., has no monodromy around 0), so that the above contour 
integral is well-defined. 

More concretely, let us write 

A('\t) = J2ji'Ht)r, ^W(^) = ^J«(^)J^ 

a a 

where {J"} is a basis of g. Let ^f{t) be the matrix elements of the fundamental solution 
of the equation (j3.8p . which may be written explicitly as follows: 

«(t) + ^j(0)(t)j;/«(t) = 0, 

c d 

where {/^^} are the structure constants of g. Then the integral in (|3.7p is equal to 

(3.9) Yl 1 4'Hmut)dt, 
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where 



b 



with the coefficients Cb{(p) determined by (p. 

Since the algebra of invariant functions on G is isomorphic to the polynomial algebra 
in ^ = rank(0) generators, the elements ^^(i), where if are invariant functions on G, 
span an ^-dimensional subspace in the space of all solutions of (13. 8p . The corresponding 
integrals (13. 7p give rise to commuting Hamiltonians of an integrable system on the space 
of connections of the form (13. 6p (in which e is taken as the spectral parameter). 

We now apply this in our case, which fits in the statement of Lemma [1] with e = z, 
= A{t) and = x- Then we obtain that the z-linear term in the expansion 

of (p{M{z)), where ip G (FunG)*^, is given by the formula 



Higher terms in the z-expansion of (p{M{z)) are given by more complicated formulas, 
involving multiple integrals of solutions of the equation 



(in the adjoint representation of g). The solutions ^{t) of this equation, and hence 
all of these Hamiltonians, are non-local, in the sense that they cannot be realized as 
integrals of differential polynomials in the matrix elements of A{t). 

Finally, observe that both local and non-local classical Hamiltonians Poisson com- 
mute with the Heisenberg-Poisson algebra generated by the Fourier coefficients of the 
Cartan components of A{t). Indeed, we obtain our Hamiltonians in effect by differ- 
entiating invariant functions on qI (namely, ip of the monodromy) in the direction of 
X € f)* CqI- Therefore these Hamiltonians are invariant under the action of the central- 
izer of x in 0) which is precisely the (homogeneous) Heisenberg subalgebra i) of g. This 
is equivalent to the fact that the Hamiltonians Poisson commute with the Heisenberg- 
Poisson algebra generated by f) (in the case of local Hamiltonians this commutativity 
is discussed from a different point of view in |FF6] ). 

We now wish to quantize these Hamiltonians. The following remark might be useful 
for understanding what follows. 

Recall that the Hamiltonians of the classical shift of argument Gaudin model for 
a finite-dimensional Lie algebra g are obtained by taking the derivatives of invariant 
polynomials on g*. Each invariant polynomial P on g* has finite degree. Therefore, 
when we shift the argument by zx and expand in powers of z, we obtain a finite 
polynomial in z, whose coefficients are the derivatives of the P in the direction of x- 
In the affine case the role of g* is played by g^, which is the space of connections (13. 3p . 
The role of invariant polynomials on g* is now played by the spectral invariants (p{M) 
of the monodromy of the connection. As functions on g^, they have "infinite degree". 
When we apply the shift of argument by zx, the nth coefficient of the z-expansion of 



(3.10) 




The Poisson commutativity ()3.5p of the Hamiltonians then implies that 
{HM^MX)} = 0, (FunG)«. 



dt^{t) + [A{t),^{t)]=0 
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the corresponding functions <p{M(z)) will be equal to the nth derivative of this function 
in the direction of x- This way we obtain the non-local classical Hamiltonians. On the 
other hand, it is instructive to think of the coefficients in the asymptotic 2:~^-expansion 
as expressions obtained by differentiating this function (oo — n) times. These are the 
local classical Hamiltonians. 

Now, the quantization of the derivatives of the invariant functions on g* in the di- 
rection of X gives rise to the quantum Hamiltonians of the shift of argument Gaudin 
model, as discussed in Section [2l As explained in [R \ IFFT| . to construct these Hamilto- 
nians we really need to use the completed enveloping algebra of q{{z)) (more precisely, 
of its critical central extension). The quantum Hamiltonians are constructed, roughly 
speaking, by taking derivatives of central elements of this algebra. Likewise, the quan- 
tization of the spectral invariants of the monodromy of the connection gives rise to the 
quantum Hamiltonians of the affine Gaudin models (quantum soliton systems). These 
Hamiltonians should also be interpreted as coming from the center of a completed en- 
veloping algebra, but now of the double affine algebra q{{z)) (or some extension of it). 
We expect that this center contains quantum analogues of the spectral invariants of 
the monodromy of the connection, and the local and non-local quantum Hamiltoni- 
ans may be obtained by taking the "derivatives" of these central elements (as in the 
finite-dimensional case), either finitely many times or infinitely many times. 

While at present we do not have a precise definition of the completed enveloping 
algebra of q{{z)) or its center, the quantum Hamiltonians of the affine Gaudin models 
do exist and may be constructed explicitly, as explained in the next section. We hope 
that by analyzing these Hamiltonians we may gain insights into the structure of the 
center of the enveloping algebra of 'g((z)). 

3.3. Quantum Hamiltonians. The quantization of the local Hamiltonians has been 
discussed in our paper |FF6| . in which we explained how to construct them using the 
free field (Wakimoto) realization of affine algebras. We refer the reader to [FF6] for 
more details. According to Conjecture 1 of [FF6j . all of the classical local Hamiltonians 
of the g-AKNS hierarchy may be quantized. Therefore we expect to have i mutually 
commuting quantum local Hamiltonians of each positive spin. 

In addition, for non-zero values of level k these quantum Hamiltonians should com- 
mute with the (homogeneous) Heisenberg subalgebra f) C g, which is the central exten- 
sion of (]((t)). 

Explicit formulas for the quantum g-AKNS Hamiltonians are unknown in general. 
But we can write down those of spin 1. First of all, we have the following formula, 
generalizing formula (|2.12p for the DMT Hamiltonians (and using the same notation): 






■+ 
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where A^_^ is the set of positive real roots of g. Here we set 



and similarly for fa{w), and assume that x is regular so that (a, x) 7^ for all a G A+. 
The normal ordering is defined as in |FBj . 

These operators mutually commute, but the drawback is that they do not commute 
with f). However, we can modify them slightly to make them commute with f) provided 
that the level k is non-zero. Namely, let us set 

(3.12) f^{x) = 

|q' y ■ (^ea{w)fa{w) + fa{w)ea{w) - ^haiw)"^^ :wdw, 7 6 f)*> 



where /iq, € is the element corresponding to the root a € f)* via the normalized 
bilinear form. 

The following lemma will be proved in the Appendix^ 
Lemma 2. 

(1) For any regular x G 1]* the Hamiltonians T^(x),7 S ^* , commute with each other. 

(2) For any regular x G f)* the Hamiltonians T^(x),7 £ f)*, commute with each other 
and with the Heisenberg suhalgehra f) C g. 

Next, we discuss the quantization of the non-local Hamiltonians. In order to do 
this, we need to construct quantum analogues of the solutions ^(t) of the differential 
equation 

(3.13) {kdt + A{t))-^{t)={), 

taking values in the adjoint representation of g. It has been shown by V. Knizhnik 
and A. Zamolodchikov in |KZj that these quantum analogues may be obtained from 
the primary fields of the WZW model corresponding to the adjoint representation. We 
will now recall how this works for an arbitrary finite-dimensional representations V of 
g. In what follows we will assume that fc is a complex number not equal to the critical 
value — /i^. 

By definition, for given highest weight modules Mi and M2 over g of the same level 
k, the primary field corresponding to is a formal power series 



•^viw) = ^ <^v,(n)W~'' : Ml y ^ M2, 



which is g-invariant, where g acts on V by evaluation at w, 



'''The operators Tj{x) may be obtained as certain limits of the connection operators studied in 
[FMTV] , and part (1) of Lemma [2] may be derived from the commutativity of these operators proved 
in [FMTV] (see Section ITSi below for more details). 



QUANTIZATION OF SOLITON SYSTEMS AND LANGLANDS DUALITY 



17 



This is equivalent to the Fourier coefficients ^v,{n) having the fohowing commutation 
relations with the affine algebra g: 

Here on the left hand side we consider the action of on Mi and M2, and on the right 
hand side the action of J" on V. 

If Ml and M2 have highest weights Ai and A2 and V = has highest weight 
jjL, then the corresponding conformal dimensions (the highest weights of the Segal- 
Sugawara Virasoro algebra) are equal to A(Ai), A(A2), and A(/i), respectively, where 
A(z^) = C{fi)/{k + h'^), and C(z^) is the scalar by which the normalized Casimir element 
of U{q) acts on V^. The true primary field is then 

As explained in |KZj . it satisfies the following differential equation: 

(3.14) {k + h'')d^^v{w) - 5^ :riw){Ja ■ ^v){wy. = 0, 

a 

which may be viewed as a quantization of the differential equation (|3.13p written in the 
representation V (note the "quantum shift" k 1— k + h^). Therefore we will consider the 
primary fields ^viw) as quantizations of solutions of the equation (|3.13p . Moreover, 
we see that the monodromy of this ^v{w) is equal to exp 27rz(A(A2) — A(Ai) — A(^)). 

Suppose now that Mi and M2 are Verma modules Mx^^k and Mx2,k- It is known that 
for generic k the space of primary fields of type V = Vf^ acting from M^a^ ^ to Mx2,k 
has the dimension equal to the dimension of the component of weight A2 — Ai in V^. 
Moreover, it is known that this space is isomorphic to the component of weight A2 — Ai 
in the irreducible representation of the quantum group Uq{Q) with highest weight fi, 
where q = exp{iTi/{k + /i^)). 

In particular, we find that if Mi = M2 = Mx, then the space of primary fields ^{w) = 
^q{w) of type Q (the adjoint representation) acting from Mx to itself is isomorphic to the 
weight subspace of q, which is the Cartan subalgebra of q. It therefore has dimension 
i, the rank of 0. It is clear that this subspace is precisely the quantization of the space 
of single-valued solutions of the equation ()3.13p . spanned by ^^p, ip € (FunG)'^. Let us 
choose a basis ^p{w),p = 1, . . . , i, oi this space of primary fields. 

Then the vector space of the quantum Hamiltonians that quantize the Hamiltonians 
Hp{x) given by formula (j3.10p is just the span of the 0th Fourier coefficients of the 
fields ^p{w), 

(3.15) Hpix) = [ ^p(-) - = / ^ 

J w J w 

of the primary fields ^p{w) : Mx — > Mx- Here we use the fact that 

A(0) = /i^/(fc + /i^), 

where 9 is the maximal root, which is the highest weight of the adjoint representation 
of 0. 
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We claim that the Hamiltonians Hp{'x),p = commute with each other and 

with the local Hamiltonians (|3.1ip . This is checked by a direct calculation which will 
be presented in [FF7J. 

The quantization of the classical non-local Hamiltonians, corresponding to the higher 
terms in the z-expansion of the monodromy ip{M{z)) is expressed by multiple integrals 
of the adjoint primary fields (^{w)w^^ /^^^^^^ (including the components acting between 
different Verma modules). We will now briefly sketch the construction of these integrals. 
The details will appear in |FF7j . 

This construction is analogous to the construction of (BLZH IBHK| . where non- 
local quantum KdV Hamiltonians were constructed using bosonic screening operators 
corresponding to the simple roots of g. These operators satisfy the g-Serre relations, 
and hence, in a certain sense, generate the quantized enveloping algebra C/g(n_), where 
n_ is the lower nilpotent Lie subalgebra of g (see [BMP! IFF5] for the precise definition 
of what this means). The non-local quantum KdV Hamiltonians are then obtained 
from elements of C/g(n_) which correspond to singular vectors in the Verma modules 
over C/g(g) of critical level. 

More precisely, as explained in |FF5j . whenever we have a singular vector of weight 
^ in the Verma module over f7g(g) of highest weight A and level 0, we have a well- 
defined operator i^x Hfj., where i:\ denotes the Fock representation of a Heisenberg 
Lie algebra with highest weight A. This operator is constructed as follows. Let us write 
the singular vector as Pv\ , where v\ € is the highest weight vector and P G Ug (n_ ) . 
This operator is obtained by substituting the bosonic screening operators 

/du) 
Si{w) — , i = 0,...,£, 
w 

where the Si{w) = e^'^'-^^^ are the bosonic primary fields corresponding to the simple 
roots in q, into P instead of the generators fi,i = 0, . . . , i. Here q = exp(7ri/3^). Now, 
if we replace the screening operator So by the operator 

(3.16) So= [ Soiw) w^'P' ^, 

J w 

and keep Si,i = 1, . . . ,i, unchanged, then the same procedure will work, except that 
we will now assign operators acting on Fock representations to singular vectors in the 
Verma modules of critical level k = —h^ instead of level k = 0. There are many more 
of those than at /c = 0, and one can check that they commute with each other. These 
are the non-local quantum Hamiltonians of [BLZl^ IBHKl IFFS] . 

The non-local Hamiltonians of the affine shift of argument Gaudin model may be 
constructed in a similar fashion. The most direct way to do it is to replace the Virasoro 
(or, more generally, W-algebra) screening operators Si with the screening operators 
for the affine algebra g arising from the free field (Wakimoto) realization of g. These 
screening operators have been defined in [FF6] . where it was shown that they satisfy 
the same Serre relations of C/q(n_). Hence, substituting them in the formulas discussed 
above, we construct non-local commuting Hamiltonians of our model. 

However, we want to stress that it is not necessary to use the free field realization 
to define these quantum non-local Hamiltonians. They may be defined intrinsically in 
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terms of the adjoint primary fields of the WZW model associated to g. In order to 
do that, we need to replace the algebra Uq{n-) by a smaller quantum algebra, which 
is a quantization of U{q-), where g_ = t~^0[t~^]. To explain this more precisely, 
let us recall that the algebra C/q(n_) may be constructed using certain rank one local 
systems on the configuration spaces (see (SV2l|Vj). Moreover, homology classes with 
coefficients in these local systems correspond to singular vectors discussed above. These 
local systems describe the monodromy properties of the above screening currents, and 
the non-local Hamiltonians may be interpreted as integrals of products of the screening 
currents over homology cycles of these local systems (see |FF5j for more details). 

Likewise, the quantization of U{g-) that we need here, which we denote by Ug{Q-), 
may be constructed using local systems on the configuration spaces that capture the 
monodromy properties of the adjoint primary fields ^g{w). These local systems are 
defined on the trivial vector bundles on the configuration spaces diag with the fibers 
and the connection given by the corresponding KZ operators. The monodromies 
of these local systems around the diagonals are given by the ii-matrices of Uq{Q) in 
the adjoint representation. Now, the commuting non-local Hamiltonians quantizing 
the coefficients in the z-expansion of (p{M{z)) are constructed by integrating products 
of the primary fields over homology classes with coefficients in these local systems, 
which in turn correspond to singular vectors in C/q(0_)-modules. Since we want to use 
singular vectors at the critical level k = —hy (and not /c = 0), we need to insert the 
factor in the integral (j3.15p . This is analogous to inserting the factor w^^^ 

in formula (j3.16p . 

The simplest of these non-local Hamiltonians are the operators Hp{x) given by for- 
mula (j3.15p . They correspond to the quantization of the z-linear term in the expansion 
of ^p{M{z)). Higher degree terms correspond to integrals over more complicated homol- 
ogy classes with coefficients in the above local system. Applying the free field realization 
of 0, we may express these Hamiltonians as contour integrals of the g screening currents 
and obtain formulas similar to those of |BLZH IBHK] . 

3.4. Bethe Ansatz. Next, we discuss the eigenvectors of the quantum Hamiltonians. 
The Bethe Ansatz described in Section [2.31 has a straightforward generalization to this 
case. The formula (|2.9p for the Bethe vector and the Bethe Ansatz equations (j2.10p 
easily generalize to the affine case. We just need to enlarge the set I to I and include 
the generator Jq (corresponding to the simple root ao of g) and view v as an element 
of f)*. Then for any collection of distinct non-zero complex numbers Wj,j = 1, . . . ,m, 
and labels ii, . . . ,im ^ I of the Dynkin diagram of g, the vector (l){w\^ , ■ ■ ■ , w^m") given 
by formula (j2.9p belongs to the Verma module over g. 

We expect that if the Bethe Ansatz equations (j2.10p are satisfied, then the vec- 
tor (p{w^i , . . . ,wl:^) is an eigenvector of the quantum Hamiltonians discussed above 
(both local and non-local). Since we expect that these Hamiltonians commute with the 
Heisenberg Lie agebra f), it is natural to assume that these vectors are also invariant 
under the Lie subalgebra tf)[[t]] C f) (we have checked that this is indeed true in the 
simplest examples). Thus, we arrive at the following conjecture. 
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Conjecture 1. // the Bethe Ansatz equations ()2.10p are satisfied, then the vector 
(p{w^i , . . . ,10^:^) is a ti}[[t]]-invariant eigenvector of all local and non-local quantum 
affine Gaudin Hamiltonians, in particular, the Hamiltonians r^,7 € \]* , given hy for- 
mula (I3.12|) . and Hp{x),p = 1, ■ ■ ■ ,i, given by formula (I3.15p . 

Note that (j){w\^ , • • • , wj^ ) is a homogeneous vector of weight A — Y^JLi (^ij ■ It is clear 
from formulas (I3.12|) and (13. lip that any homogeneous tf)[[t]] -invariant eigenvector of 
r-Y,7 € f)*, is automatically an eigenvector of T^,7 € f)*. Therefore Conjecture[T]implies 
the following statement. 

Conjecture 2. The vector (j){w^^ , ■ ■ ■ , w^m) eigenvector of the affine DMT Hamil- 

tonians (13. lip if the Bethe Ansatz equations (I2.10p are satisfied. 

To prove this conjecture, we cannot use the argument used in [FFTj because we do 
not know that the affine DMT Hamiltonians arise by the same mechanism as the DMT 
Hamiltonians in the finite-dimensional case (that is, from the center of an enveloping 
algebra of the affinized Lie algebra). But this conjecture can probably be derived from 
the results of |FMTV] (see Section [4.81 below for more details). 

The next question is to describe the common eigenvalues of the quantum Hamilto- 
nians that occur on highest weight g-modules, such as Verma modules and irreducible 
integrable representations. In particular, we wish to describe the eigenvalues on the 
Bethe eigenvectors. 

According to Theorem [2l in the finite-dimensional case the common eigenvalues of 
the quantum shift of argument algebra are encoded by ^g-opers on satisfying 
certain conditions. In particular, the common eigenvalues on the Bethe vectors (|2.9p 
are encoded by the ^g-opers obtained by applying the Miura transformation to the 
connection (12. lip . We wish to generalize these results to the affine case. This will be 
done in the next section. 

4. Affine opers 

The notion of g-opers, where g is a finite-dimensional simple Lie algebra, was in- 
troduced by Beilinson and Drinfeld [B Dll IBD2] following Drinfeld and Sokolov |DS] . 
Denote by G the Lie group of inner automorphisms of the Lie algebra g and B its Borel 
subgroup. A Q-oper on a smooth curve X is by definition a triple {3',V ,9'b), where 
3' is a principal G-bundle on X, V is a connection on 9" and 3'b is a i3-reduction 
of 3' such that locally on X, if we choose a local coordinate z and trivialize ^b, the 
connection acquires the form 

e 

(4.1) V = a, + ^^,(z)/, + v(z), 

i=l 

where each ipi{z) is a nowhere vanishing function, and v{z) is a b-valued function. 

In this section we generalize the notion of g-opers and define g-opers, where g is 
an affine Kac-Moody algebra, following iF2]. We then use these affine opers to give a 
conjectural description of the spectra of the quantum Hamiltonians of the affine Gaudin 
model introduced in the previous section. 
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4.1. Definition. We now consider another copy of the affine Kac-Moody algebra g, 
which we will use to describe the spectra of the quantum affine Gaudin Hamiltonians. 
It is important to realize that this affine Kac-Moody algebra has nothing to do with 
the affine Kac-Moody algebra used in the construction of the Gaudin models above 
(in fact, we expect that the two are Langlands dual to each other). To emphasize this 
difference, we will now use A as the loop variable. Thus, in this section we set 

Q = Cd\xg, 

where 

^ CI ^ ^ 0((A)) ^ 0. 

It will be convenient to use the following slightly unusual convention for the Cartan 
decomposition: 

= n_ e^en+, 

where 

n+ = (n+ 1) e (0 A-^C[A-^]), 
n_ = (n_ 1) e (0 AC[[A]]), 

^= ([) (g) 1) e CI e Cd. 

Now let G be the Lie group (more precisely, ind-group scheme) associated to (see, 
e.g., [Ka| for the precise definition). We will denote by G the quotient of G by its 
center. Thus, the Lie algebra of G is 

1 = Cdx0((A)). 

The group G comes with the lower unipotent and Borel subgroups and B- (which 
are proalgebraic groups) corresponding to the Lie subalgebras 

n_ = n_ and b_ = (f) Cd) n_ 

of 0, respectively, and the upper unipotent and Borel subgroups and (which 
are group ind-schemes) corresponding to 

n+ = n+ and b+ = ({) Cd) n+, 

respectively. 

Let X be a smooth algebraic curve. We will use the notion of principal G-bundle 
on X and a connection on it given in |F2j . Roughly speaking, such a bundle may be 
trivialized on open subsets of X (in either Zariski or analytic topology). Having chosen 
such a trivialization and a local coordinate z on a an open subset U, a connection is 
by definition a first order differential operator V = c^^ + A{z), where A{z) G 0. If we 
change the trivialization by g{z) G G, V transforms by the usual formula 

V ^ + gA{z)g-' - {d,g)g-\ 

Under a change of coordinates z = f{s) the operator V also transforms in the usual 
way: 

V^ds + ip'{s)A{ip{s)). 
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Now, by analogy with the finite-dimensional case [DS^ IBDlj IBD2j . a Q-oper on X 
may be defined as a triple (3", V, S"-^^), where 3" is a principal G-bundle 3' on X, V is 

a connection on 3' and 3"-^^ is a -reduction of 3", such that locally, with any choice 
of local coordinate z and trivialization of 3"-^^ , the connection operator has the form 

e 

(4.2) ^ = d^ + Y,Mz)fi + ^{z), 

i=l 

where each ipi{z) is a nowhere vanishing function, and v(z) is a b+-valued function. 

We note that another version of affine opers, introduced earlier in |BFj . has b+ and 
b- switched. For comparison of the two definitions and the explanation as to why the 
definition given here is the one relevant to the Gaudin model, see |F2) . 

However, this naive definition appears to be inadequate. The reason is that for ind- 
groups, such as G, a connection does not necessarily give rise to a trivialization of 3", 
even locally analytically. This is because the equation 

{d, + A{z))^{z)=0, 

where A{z) takes values in g ~ Cd x 0((A)), does not a priori have local solutions with 
values in G ~ x G((A)) (we may well obtain power series in A unbounded from below, 
which are not well-defined as elements of G). This is related to the fact that we do not 
have an exponential map from the Lie algebra g to the group G. This means that the 
usual correspondence between fiat connections on G-bundles and G-local systems does 
not exist. In particular, the notion of monodromy of an abstract oper as defined above 
does not make sense. 

Our experience with opers suggests that their monodromy properties are important in 
the context of Gaudin models (see |F2tlF3| ). In order to have well-defined monodromy, 
we will assume that, locally on X, the oper bundle 3" admits a reduction to a proalgebraic 
subgroup of G, such as G[[A]], which is preserved by the connection V. Then the 
connection V does give rise to local analytic trivializations of this G[[A]]-subbundle, 
and hence of the G-bundle 3" as well, so the monodromy of the oper connection is well- 
defined. Thus, from now on by a g-oper on X we will understand a triple (3", V, J'-g ) as 
above such that, locally on X, 3" admits a horizontal reduction to G[[A]] or to another 
"maximal compact" subgroup of G((A)). We denote the set of g-opers on X by Opg(X). 

Let 

P-i = ^ ft = ^ fi + ee ® A. 

We will say that a g-oper on X\x has regular singularity at x if the connection operator 
may be brought to the form 

(4.3) V = a, + i(p_i+v(z)), v(z) G b+[[z]], 

on the punctured disc near x, were z is a local coordinate at x. We define the 
residue of such an oper as in the finite-dimensional case (see, e.g., |F2j ). so that it takes 
values in f)/Waff- We will denote by the projection f) — > f)/VFaff. 
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Let u he a dominant integral coweight such that Vi = {v, 6ti) G Z+ for all i ^ I. We 
will say that a g-oper on X is u-regular at a; G X if it may be brought to the form 

(4.4) v = a, + ^z'^»/, + v(z), 

where v(z) G fi+il^;]], on . 

We also define opers with singularity of order 2 and their 2-residue similarly to the 
finite-dimensional case (see [FFT] ). 

4.2. Miura opers and Miura transformation. Let H be the subgroup of G corre- 
sponding to the Lie algebra = Cd© [). Let p be the element of f) such that {oi, p) = 1 
for all i £ I. It defines a homomorphism — > if. Let be the -ff-bundle on X 
obtained as the push- forward of the -bundle on X corresponding to the canonical 
line bundle on X. Connections on fi'' are described as follows. If we choose a local 
coordinate z on an open subset of X, then we trivialize fi'' and represent the connection 
as an operator 

(4.5) V = d^ + u{z), u{z)eij. 

If s is another coordinate such that z = (p{s), then this connection will be represented 
by the operator 

(4.6) ds + ^'{s)u{ip{s))-p-^. 

We will call them Cartan connections. 

We define the Miura transformation from Cartan connections to g-opers taking an 
operator ()4.5p to the gauge equivalence class of 

(4.7) V = d,+p^i + vi{z). 

The Miura transformation may be understood more conceptually in terms of the 
Miura opers. By definition [F4], a Miura g-oper on X is a quadruple (9", V, ^''g^ 1 ) > 
where (9", V, 9"-^^) is a g-oper on X and is a i3_-reduction of 3" which is preserved 
by V. It is called generic if the reductions "3^^^ and are in generic relative position 
(see |F4j ) . We have a map from Cartan connections to generic Miura opers taking an 
operator ()4.5p to the quadruple in which 3" is the trivial bundle with the connection 
(j4.7p and the tautological reductions '3^^j^-,'3^^_ satisfying the above conditions. This 
map sets up a bijection between generic Miura opers and Cartan connections (see 
[F4j ) . The inverse map sends the quadruple as above to the -ff-bundle 3"-^ /N- with 
the connection induced by V. 

An arbitrary Miura oper on X is always generic away from finitely many points of 
X. At each of the remaining points the reductions 3"-|j^ and SF^ are in relative position 

labeled by an element tt; 7^ 1 of the (affine) Weyl group W^s of G. The corresponding 
Cartan connection has a singularity at such a point with the residue equal to p — w{p) 
(see |F2] for more details). 
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The Miura transformation is thus defined as the composition of the map from the 
Cartan connections to generic Miura opers and the forgetful map from Miura opers to 
opers. 

4.3. Parabolic Miura opers. Finally, wc introduce the notion of parabolic Miura 
oper. By definition, a parabolic Miura g-oper on X is a quadruple V, 3^;^^ , ), 

where (if, V, -Fg^) is a g-oper on X and 5^ is a reduction of 5" to G_ = x G[[A]] 
which is preserved by V. We will call it generic if S"^^ and S'q_ are in generic relative 
position. We then have the following map from the space Opg(X) of g-opers on X 
to the space of generic parabolic Miura g-opers: it sends a g-oper represented by the 
connection 

V = a, + J]/i + v(z), v(z)Gb+, 

to the quadruple in which is the trivial bundle, 

V = V + /o, 

and -Fq are the tautological reductions of 9". However, this map is not bijective, 

because we are missing the one-dimensional subgroup of G- corresponding to the 
element d G [}. 

In order to obtain a bijection, consider the Cartesian product of Opg(X) and the 

space Connx(^^''), where h = {S, p) is the Coxeter number of g. Define a map from 
Opg(X) X Connx(^^'^) to the space of generic parabolic Miura g-opers sending a pair 
consisting of an oper V € Opg(X) and a connection v' G Connx(f^'*) of the form 

to the connection 

V = V + /o + 7(^K 
This map is already bijective. To construct the inverse map, observe that 

(where G^^^ is the first congruence subgroup of G[[A]]). To a quadruple (?", V, 9^^^ , 3^q_ ) 

as above, we associate a G x C^-bundle "J/G^^^ together with a connection and a B- 
reduction% /G^'^\ This bundle is a product of the g-oper bundle and the -bundle 

corresponding to Therefore these data give rise to a point in Opg(X) x Connx(^^'*)- 
Taking the composition with the forgetful map from parabolic Miura opers to opers, 
we now obtain a map from Opg(X) x Connx(^^'*) to g-opers, which is a parabolic version 
of the Miura transformation. This gives us a convenient way to represent g-opers in 
terms of g-opers, which we will exploit in what follows. 

Note that any parabolic Miura g-oper on a curve X is generic away from finitely 
many points of X. At each of those points the relative position of ffg^ and 3^q_ is 
measured by a coset w G Ty\Waif , and the corresponding g-oper has regular singularity 
with the residue w{—w{p)), where the bar stands for the projection f) ^ f) along Cd 
and vj is the projection [) — >■ \)/W^ as before (note that w{—wi{p)) = w{—W2{p)) if 
w\ ~ W2 in VF\VFair, so this formula is well-defined). 
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Remark 1. There is an analogous notion of parabohc Miura oper in the finite-dimen- 
sional case. Let P be a parabolic subgroup of G. Then we define a Miura (g, P)-oper 
to be a quadruple (3^, V, IFb, ff'p), where {3',V,3^b) is a g-oper (see the beginning of 
this section) and S'p is a P-reduction of 3" which is horizontal with respect to V. (In 
the affine case the role of P is played by We have the obvious notion of generic 

Miura (g, P)-oper. 

Let L = P/[P,P]. Then we define, in the same way as above, a bijective map from 
the space of L-opers to generic Miura (g,P)-opers. Composing it with the forgetful 
map from Miura (g, P)-opers to g-opers, we obtain a parabolic analogue of the Miura 
transformation map from L-opers to g-opers. The ordinary Miura opers correspond to 
the case P = B. □ 

4.4. Spectra of quantum Gaudin Hamiltonians and afRne opers. Now let ^g 

be the Langlands dual Lie algebra of g, whose Cartan matrix is the transpose of the 
Cartan matrix of g (thus, ^g is a twisted affine Lie algebra if g is non-simp ly laced). 
We will use the upper index L to denote all of the above groups corresponding to ^g, 
such as ^f), ^H, etc. In particular, we have a canonical identification ^f) = and so 
we will consider a weight z/ € f)* as an element oi ^\). By abusing notation, we will 
denote by the same symbol its projection onto ■^f) = ^I)/C1 = Cd. 

As in the finite-dimensional case (see [FFT] and Section ES]), we attach to each 
solution of the Bethe Ansatz equation ()2.10p . where ij € /, the Cartan ^iL-connection 

m 

(4.8) v = a,-x-- + E^^- 

Z ^ Z — Wj 

In the same way as in [F2], Lemma 2.10, we prove the following: 

Lemma 3. The Bethe Ansatz equations (|2.10p are equivalent to the condition that the 
^Q-oper V = V + p_i obtained by applying the Miura transformation to the connection 
(j4.8p is regular at the points Wj,j = l,...,m, or, equivalently, has no monodromies 
around these points. 

Note that this ^g-oper has regular singularity at the point z with the residue and 
cxD with the residue determined by i^, and order 2 singularity at oo with the 2-residue 

■^ix)- 

This suggests that the true parameters for the common eigenvalues of the quantum 
Gaudin Hamiltonians are not the Cartan connections (I6.15p . but the corresponding ^g- 
opers obtained by applying the Miura transformation. Thus, we propose the following 
description of the spectra of the quantum affine Gaudin Hamiltonians (including the 
affine DMT Hamiltonians T^{x) and the non-local Hamiltonians Hp{x))- 

The common eigenvalues of the quantum affine Gaudin Hamiltonians on the 
subspace of t\][\t\\-invariant vectors in an irreducible Q-module with highest 
weight v and level k ^ over g are encoded by ^q- opers on with regular 
singularity at with residue w{—v — p) and singularity of order 2 at oo with 
2-residue 7r(x). 
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Here it is important to note that the projective hne P^, on which these ^g-opers are 
defined, has a global coordinate z, which is nothing but the spectral parameter of the 
affine Gaudin model (in particular, z is the spectral parameter in the L-operator p.4p ). 

4.5. The case of sl2. What do these ■^g-opers look like? Consider the case when 
g = ^0 = 5(2. The Cartan subalgebra ^i) has a basis consisting of 



h 



1 
-1 



1, and d, and ^f) is spanned by h and d (recall that d acts as Xdx). Without loss of 
generality we may, and will, consider only highest weights u such that {u, d) = 0. They 
are determined by the numbers k = {i^,!} and £ = ^{v,h). Then the corresponding 
element oi^\) = Ch © Cd is represented by ih + kd. 

We have two simple roots, ai and ao, which in^t) are represented as ai = h and ao = 
—h. We therefore divide the set of points {wi, . . . , Wm} into two subsets: {wl , . . . , } 
and {w^, . . . , tf^oi- "^^^ connection (j4.8|) gives rise to the Miura sl2-oper 

(4.) a..(t' 

where 



mi ^ mo ^ 



(4.10) u{z) = -x-'- + Y^ r - E ■ 

z ^-^ z — w ■ ^-^ z — w - 

j=i J j=i 3 

Concretely, the statement of Lemma [3] is realized as follows. The connection in (j4.9p 
has apparent singularities at the points z = Wj and z = w'j (in addition to z = 
and z = 00), but the condition of Lemma [3] is that these singularities are "fake"; 
that is, they can be removed by a gauge transformation (by an element of -B+). This 
imposes certain conditions on the points wj and w^, which are in fact the Bethe Ansatz 
equations ()2.10p . 

Let us see more precisely how this works. This is in fact a local question, which 
involves only the restrictions of the Miura oper (j4.9p to the formal discs around these 
points. 

Suppose that we have a Miura 5[2-oper of the form 

(4.11) 4 +("</) _4,)+x(.)<i 
on the formal disc around a point w such that 

u{z) = ^ Un{z - w)"-, x{z) = ^ Xn{z - wj" . 

n>-l n>0 

Suppose first that u_i = 1. This is the case when u{z) is given by formula (j4.10p and 
w = wY Let us apply the gauge transformation by 

i_ 

.0 1 



QUANTIZATION OF SOLITON SYSTEMS AND LANGLANDS DUALITY 



27 



Then we obtain the operator 

2 



5.+ 



w 



u{z) 



w 



1 




+ 



where the dots stand for the terms regular at z = w. Therefore this operator is regular 
at z = w a and only if the equation 



(4.12) 



-uo = 



is satisfied. 

Next, consider the case when n_i = 
()4.10p . Applying the gauge transformation by 

1 



1, which corresponds to w = Wj in formula 



we obtain the operator whose expansion near z = w has the form 



9.+ 



w 



-u{z) 



1 



x{z) 



z — w 







+ ..., 



where the dots stand for terms regular at z = w. Therefore this operator is regular at 
z = w ii and only if the equation 



(4.13) 



Uo + 







is satisfied. 

Thus, we obtain that the connection (|4.9p with u_i = 1 (resp., u_i = —1) is gauge 
equivalent to a connection regular at z = it; if the equation (j4.12p (resp., (|4.13p ) is 
satisfied. Conversely, one shows, following the argument of [F2j . Lemma 2.10, that 
if this equation is not satisfied, then the connection (j4.9p has non-trivial monodromy 
around w. Therefore it is impossible to make it regular by a gauge transformation. 



Let us apply this in the case of our Miura oper (j4.9p . Then equations (|4.12p for 
w = Wj become 



(4.14) T + Y.— 0+^ = 0' j = l,...,rni, 

W!: ^ W\ — W% — 
J s^] J * s = l J * 

and equations (j4.13p for w = Wj become 

k _ £ rni ^ 1 

(4.15) ^^ + Y.— I-Y.— 0-^ = 0' j = !,■■■, mo. 

3 s=l J * s^j 1 

These are precisely the Bethe Ansatz equations ()2.10p . 
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4.6. Second order operator. Now we wish to rewrite (14. 9 p as a second order differ- 
ential operator. To this end, we apply the gauge transformation by 



(4.16) 



1 -u{z) 
1 



Then we obtain the following operator: 

(4.17) 4+(; "(y^)4<i. 

where 

v{z) = u{z)'^ + dzu{z) 

, . 2 i . ^k\\ £(^ + 1) ^2 ^ A: 1 

I wi \ z z^ (z — w^Y vr- z — w^- 

\ j=i 3 j j=i ^ j=i -J ^ 

In deriving this formula, we have used the Bethe Ansatz equations (|4.14p ~ (j4.15p . They 
imply in particular that vi^z) has no singularities at tnj, j = 1, . . . ,mi. 

In addition, the Bethe Ansatz equations ()4.15p imply that the coefficients Vj.^ in the 
expansion of vi^z) at z = Wj, 



2 

(z - w^y 

^ J' n>— 1 

satisfy the algebraic equation 

\ ( k 

(4.19) 7 —0 o^i.o + = 0' J = !,•••, "io- 

\ j ) j 

Note also that Vj-i = k/w^. 

For instance, if A: = 0, then we have vj-i = and the equation (|4.19p reduces to 
Vj^i = for all j = 1, . . . ,mQ. Hence we obtain that in this case the condition is that 
the coefficients in front of {z — Wj)~^ and {z — w^) in the expansion of v{z) in Laurent 
series in {z — w'j) should vanish for all j = 1, . . . , mo. 

These equations guarantee that on the formal neighborhood of each of the points 
w^,j = 1, . . . ,mo, the s[2-oper (j4.17p . viewed as a connection on the trivial G-bundle 
on the punctured disc around this point, is gauge equivalent to a connection without 
singularity. In other words, this 512-oper has no monodromy around the points Wj,j = 
l,...,mo. 

This is in fact a local condition on this oper which involves only its restrictions to 
the formal discs around these points. Let us analyze this local condition in more detail. 

Suppose that we have an 5(2-oper of the form 

(4.20) ^^'\-'')+xiz)d, 
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on the formal disc around a point w such that 

(4.21) V{z) = . ^ ,0 + Vn{z-wr, x{z) = y^Xniz-wr. 

^ ' n>-l n>0 

What is the condition that the connection ()4.17p on the trivial G-bundle on the punc- 
tured disc around the point w (where G = >< G((A))) is gauge equivalent to a 
connection without singularity? 

This connection preserves the trivial G_-bundle, where G- = xG[[A]]. Taking the 
quotient modulo the subgroup k where G^^^ is the first congruence subgroup 
of G'[[A]], we obtain the connection 



(4.22) + 



v{z) 

1 



on the trivial G-bundle, which we view as an sl2-oper. If the connection (I4.20p has 
no monodromy around z = u), then the same is true for the connection (I4.22p . This 
implies that the following equation on the coefficients in the expansion of v(z) at z = w 
must hold: 

(4.23) - vov-i + f i = 

(see, e.g., |FFTj . Sect. 4.4). It is easy to check that this equation is satisfied if and 
only if v{z) may be represented as the Miura transformation 

(4.24) v{z) =u{zf -d,u{z), 
where 

U{z) = + V Un{z - W)"". 

z — w ^-^ 

n>0 

Thus, our connection (|4.20p is gauge equivalent, under the gauge transformation 
(j4.16p . to a connection of the form ()4.1ip with u-i = —1. As explained in the previous 
section, the latter is gauge equivalent to a connection without singularity aX. z = w \i 
and only if equation (|4.13p is satisfied. Substituting this into formula ()4.24p . we find 
that 

(4.25) u_i = xq. 
As the result, we obtain the equation 

(4.26) -xi - VQXQ + vi = 0. 

Thus, we obtain that the condition for the operator (j4.20p to be monodromy free is 
given by two equations: (|4.26p and (14.25[) . Applying this to the s[2-oper (j4.17p on 
with w = w^, we obtain the equations (|4.19p and vj-i = k/wj. 

Recall that we obtained equation (I4.19P originally under the assumption that v{z) 
may be represented globally, on the entire P^, as the Miura transformation of (14. 9p . But 
now we see that we do not need to make this assumption. In fact, the local condition 
of triviality of the monodromy around a point z = w of any 5[2-oper ()4.20p with v{z) 
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and x{z) given by formula (I4.2ip is expressed by the equations ()4.23p and ()4.25|) . which 
give rise to the equation (I4.19[) in the case of the oper (|4.17p . 

Applying the (formal) gauge transformation by z^'^ to (j4.17p , we obtain the operator 

^^+(1 

which may be rewritten as the following second order operator with spectral parameter 

mo 




.^ .2 2 „ ^ k \1 ^2 ^ A; 1 
4.27 2^x + E— oW-E— 0-^^ 



k 

— AZ 



Now we may treat A as a complex parameter. The equations (j4.19p then imply that the 
operator (j4.27p has no monodromy around the points = 1, . . . ,mo, for all values 
of A. 

Let us summarize: a connection dz + u{z) given by formula (j6.16p with the numbers 
w^jiw] satisfying the Bethe Ansatz equations (|4.14p - (j4.15p gives rise to an s[2-oper 
(j4.17p on (equivalently, the second order operator (j4.27p ) such that the coefficients 
in its expansion at w^^,] = 1, . . . ,mo, satisfy the equations (j4.19p . As we have seen 
above, these equations are the "no monodromy" conditions which may be viewed as 
analogues of the Bethe Ansatz equations (j4.15p . Our proposal for = is that for 
generic x and v the common eigenvalues of the quantum affine Gaudin Hamiltonians 
on the Verma module over s[2 are encoded by 5(2-opers of this form. 

A solution of the Bethe Ansatz equations (I4.14l )- (j4.15p gives us a solution of the 
equations (|4.19p . but we do not expect the converse to be true in general unless we 
impose additional conditions on v{z). For instance, suppose that v[z) is of the form 
(j4.18p and that in addition to having no monodromy around the points w^, which is 
ensured by the equations ()4.19p . the operator ()4.17p has monodromy ±1 around z = 
(this means, in particular, that £ is a non- negative half- integer) . Then we believe that 
such v{z) does come from a solution of the Bethe Ansatz equations (I4.14|) - (l4.15p . This 
is a conjecture, which is closely related to Conjecture [3] below. 

Note that the operator (|4.17p also gives rise to a parabolic Miura s[2-oper corre- 
sponding to the 5[2-oper 

'0 v{zy 
1 



dz + 



and the connection dz + \ 'u\ Connpi(r22), in the sense of the definition of Section [4.31 
The singular points w^j,j = 1, . . . ,mo, are the points where the two reductions, liF-^^ 
and of this parabolic Miura oper are in the relative position sq ^ W^V^^affj where 
So € W^aff is the simple reflection corresponding to the simple root uq. This is the 
generic situation. For special values of v and x tl^s positions of some of the poles 
may coalesce and these opers may degenerate further. In other words, there may 
appear points on at which the relative positions of the two reductions are given by 
elements of VF\Wafr other than sq. If the relative position at a point w is given by the 
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element so('SiSo)" G W\Wa.s,n € Z, then the expansion of the function v{z) appearing 
in (j4.18p around w will have the form 

, , (2n + l)(2n + 2) 

[z — wy 

The coefficients in the expansion of this function at w will also have to satisfy an 
analogue of equation ()4.19p . which means that our 5(2-oper has no monodromy at w. 

In addition, some of the points may tend to 0. This happens when the Verma 
module My becomes reducible. The common eigenvalues encoded be these degener- 
ations of the opers ()4.17p should correspond to the eigenvectors which appear in a 
submodule of My. 

Finally, note that we may rewrite (j4.17p as the following partial differential operator: 

d, + -\dA -v{z)-X. 

This formula underscores the fact that we are dealing here with two independent vari- 
ables: z and A. In our analysis we have treated z as a variable, but A as the parameter 
of the loop algebra 0((A)). Thus, we have considered an ordinary differential operator 
(j4.17p , in z only, but with values in an infinite-dimensional Lie algebra; namely, the Lie 
algebra Cd 0((A)). However, considering it as a partial differential operator, in both 
z and A, might also be useful in some applications. 

4.7. Integrable representations. Now we consider the common eigenvalues of the 
affine Gaudin Hamiltonians on an integrable representation of g. Recall that for each 
integral dominant weight z^, i.e., such that (i^, Oj) G Z+,i G I, so in particular k = 
{ly, 1) G we have an irreducible integrable representation V^, which is the irreducible 
quotient of the Verma module M^. By analogy with the finite-dimensional case (see 
|FFT| ■ Theorem[2l(2) and the conjecture after it), we propose the following description 
of the common eigenvalues of the affine Gaudin Hamiltonians on the space of 
invariant vectors in integrable representations of g. Note that this space is a module 
over the corresponding generalized parafermion algebra, which appears in the space of 
states of the (g, f}) coset model. 

Conjecture 3. The set of common eigenvalues of the affine Gaudin Hamiltonians on 
the space of ti)[[t]]-invariant vectors in an integrable module Vy overg (counted without 
multiplicities) is in bijection with the set of ^q- opers on ¥^ which are u-regular at 0, 
and have singularity of order 2 at oo with the 2-residue 7r(x). 

Note that all of these opers necessarily have trivial monodromy on . 



If the projection of the Bethe vector (j6.3p onto is non-zero, then the corresponding 
common eigenvalues of the affine Gaudin Hamiltonians are recorded by the ^g-oper 
obtained by applying the Miura transformation to the Cartan connection ()6.15p . It 
then satisfies the conditions listed in Conjecture [3l This provides the ffist positive test 
for this conjecture. 

It is natural to ask whether the Bethe Ansatz is complete in this model for generic 
values of x- As in the finite-dimensional case (see [FFT] ). we expect that it is complete. 
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which gives us hope that the Bethe eigenvectors may furnish an exphcit family of bases 
for all integrable representations of affine Kac-Moody algebras (parameterized by x)- 
If the Bethe Ansatz is incomplete, then there are eigenvectors other than the Bethe 
vectors introduced above. In this case we conjecture that the common eigenvalues of 
the affine Gaudin Hamiltonians would still be encoded by the ^J-opers as in Conjec- 
ture [31 However, the eigenvalues corresponding to the non-Bethe eigenvectors would 
be encoded by opers that cannot be expressed as the Miura transformation of a con- 
nection of the form (j4.8p . The structure of these opers in the finite-dimensional case is 
explained in |F3| . Sect. 5.5. We expect that the same pattern holds in the affine case 
as well. In particular, we expect that the opers corresponding to the non-Bethe eigen- 
vectors may still be expressed as the Miura transformation, but of Cartan connections 
of the more general form 

(4.28) a.-x + ^^^^^ + E^^^' 

Z ^-^ Z — Wj 

■j=l J 

where p G ^f) is defined as in Section and y and yj,j = 1, . . . , m, could be arbitrary 
elements of the affine Weyl group corresponding to ^g. They measure the relative 
position between the two reductions of the oper bundle, the oper reduction 3^b+ and 
the horizontal reduction 3'b- at the points z = and z = wj, respectively. 

Note that if y = Sj is the simple reflection corresponding to the ith simple root, 
then p — y{p) = ai. Thus, if y = 1 and all the y^-'s are simple reflections, then (j4.28p 
becomes (j4.8p . We expect that this is the generic situation. However, for special 
values of parameters some of the common eigenvalues may be encoded by the opers 
corresponding the Miura transformation of the more general connections (j4.28p . For 
instance, for g = 5(2 the weights p — yj{p) have the form v = 2nj, where nj is an integer. 
The corresponding operator (j4.27p will then have the leading term —nj{nj + l)/{t — Wj)'^ 
(formula (I4.27P corresponds to the case when all rij = 1). 



4.8. Generalization to other Kac Moody algebras. It is natural to ask which 
of the structures discussed above may be generalized to an arbitrary symmetrizable 
Kac-Moody algebra g. Let us fix x £ f)* such that 

(a,x)y^O, VaGA+, 

where A_|_ is the set of positive roots of g. The following analogues of the DMT 
Hamiltonians (I2.12p are well-defined on any highest weight module over g: 

(4.29) r,(x)= r^/i'^ei^ 7Gr, 

QGA+ 

where {ca^} and {fa^} form bases of the root subspaces go and Q-a, respectively, which 
are dual with respect to a non-degenerate invariant inner product on g. These operators 
coincide with the connection operators of a connection studied in [FMTV] . The flatness 
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of this connection, established in [FMTVj for an arbitrary symmetrizable Kac-Moody 
algebra, implies that these operators mutually commute for different 7 (and fixed x)lj 

However, note that for the affine Kac-Moody algebras these Hamiltonians are differ- 
ent from the Hamiltonians (j3.1ip . While in formula ()4.29p the summation is over the 
set of all positive roots of g, in formula (j3.1ip it is over the set of positive real roots (and 
X, 7 are elements of the dual space to the Cartan subalgebra of the finite-dimensional 
Lie algebra g, not the extended Cartan subalgebra of g). The operators (j3.1ip may be 
obtained from (j4.29p if we choose 7 to be orthogonal to the imaginary root 6 and then 
take the limit when x becomes orthogonal to 5 as well. 

Consider the Bethe vectors (f){wY , . . . ,wl^) given by formula (|2.9p . where ij G /, the 
set of vertices of the Dynkin diagram of g. We then have an analogue of the statement of 
Conjecture [2] for the operators (|4.29p . It may be derived from the results of |FMTV| by 
considering the asymptotics of the solutions of the dynamical KZ equations constructed 
in [FMTVj near the critical level, along the lines of |RV] ^ By taking the limit described 
in the previous paragraph, one can then prove Conjecture [2l 

On the other hand, as explained in }F2] . the notions of opers, Cartan connections 
and Miura transformation generalize to arbitrary symmetrizable Kac-Moody algebras 
in a straightforward fashion. Moreover, to each solution of the Bethe Ansatz equations 
()2.10p we may attach a Cartan connection by formula (14. Sp . Then its Miura transfor- 
mation is an ^g-oper which satisfies the same properties as in the finite-dimensional 
and affine cases. Therefore it is tempting to assume that analogues of our conjectures 
relating ^g-opers to the spectra of commuting Hamiltonians (in particular, an analogue 
of Conjecture [3]) also hold in this general context. 

However, for Kac-Moody algebras other than finite-dimensional or affine, it is not 
clear how to construct other commuting quantum Hamiltonians. So at present it is not 
clear what is the commutative algebra whose spectra are described by these ^g-opers. 



5. Quantum KdV system 

In this section we consider another example of the shift of argument affine Gaudin 
model, which corresponds to the g-KdV hierarchies. In this case the shift x is the 
nilpotent element of g C g corresponding to the maximal root. The corresponding shift 
of argument algebra is not so interesting by itself, so we look instead at its image in 
the Hamiltonian reduction with respect to the centralizer of X; which is n((t)). Clas- 
sically, this Hamiltonian reduction has been introduced by Drinfeld and Sokolov [DSj . 
and the corresponding commuting local classical Hamiltonians are the Hamiltonians of 
the g-KdV hierarchy, as explained in [ DSj . These Hamiltonians generate a commuta- 
tive subalgebra in the Poisson algebra obtained by the Drinfeld-Sokolov Hamiltonian 
reduction from the algebra of functions on g*. This Poisson algebra is known as the 
classical W-algebra. Its quantization is the quantum "W-algebra 'W(g) introduced in [FLj 
for g = sin and in [FF3j for general g (see also |FF2l IFF5t IFBj ). Thus, the quantum 



®This commutativity may probably also be proved by generalizing the argument used in [TL] . We 
thank V. Toledano Laredo for a useful discussion of this point. 
'^We thank A. Varchenko for pointing this out. 
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Gaudin Hamiltonians in this case are commuting elements of W(0). They are known 
as the local quantum KdV Hamiltonians |FF5] . 

As in the case of regular semi-simple shift of argument, discussed in the previous 
section, there are actually two types of Hamiltonians: local and non-local. They are 
obtained by expanding invariant functions of the monodromy at z = oo and z = 
0, respectively. The existence of the local quantum Hamiltonians was established in 
|FF4l IFF5] ■ though explicit formulas for them are unknown in general. The non-local 
quantum Hamiltonians have been constructed in [BLZ1| . Their spectra on irreducible 
highest weight modules over the Virasoro algebra have been studied in the series of 
papers |BLZ1] - [BLZ5] . In particular, in [DTH [BLZ4t IBLZ5] a remarkable conjecture 
was made relating these spectra to certain Schrodinger operators. In this section we 
show that this conjecture fits in perfectly as a special case of our general proposal that 
the spectra of quantum Hamiltonians in the affine Gaudin model associated to g are 
given by opers for the Langlands dual affine Lie algebra ^g. 

5.1. The classical model. The phase space of the classical model is defined in the 
same way as in Section 13.21 It consists of the operator^ 

(5.1) L = dt + A{t) + xz, A{t) G 0((i)) 
where now 

X = eg C 11+ C C 0((t)) 

is a nilpotent element corresponding to the maximal root of 0. The centralizer of x is 
therefore the Lie algebra n+ ((t)) , As explained in Section 12.11 instead of constructing 
a Poisson commutative subalgebra in the Poisson algebra of functions on (j5.ip . we 
may construct a commutative subalgebra in its Hamiltonian reduction with respect to 
n+((t)) and its character. We choose the character ip which corresponds to the element 

e 

(5.2) p_^ = ^f, e n_ C0((t)). 

1=1 

The resulting Hamiltonian reduction is the Drinfeld-Sokolov reduction |DSj . 

The reduced phase space of the Drinfeld-Sokolov reduction is therefore the quotient 
of the space M(0) of operators of the form 

(5.3) L' = dt+p-i + v{t), v(i) e b+((t)), 

by the gauge action of the group N^{{t)). Here 

e 



In what follows, in order to simplify the exposition, we will assume that is untwisted. But our 
discussion applies to twisted affine algebras as well. 
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is a "principal nilpotent element" of q{{z)). According to |DS| . the action of A^+((t)) on 
M(g) is free. The resulting quotient spac^EI 

M(0)=M(0)/iV+((t)) 

is isomorphic to the space of g-opers on the punctured disc = SpecC((t)). The 
Poisson algebra of local functionals on M,(q) is known as the classical W-algebra 'W(g). 
For example, M(sl2) is the quotient of operators of the form 

Q^^l^ait) m+^^^y a{t),b{t)eC{{t)), 

by the upper triangular gauge transformations depending on t (but not on z). It is 
easy to see that each orbit contains a unique operator of the form 

and hence we may identify M(s[2) with the space of such operators, or, equivalently, 
with the space of second order differential operators with spectral parameter 

d^-v{t)-z, v{t)€C(it)). 

Similarly, the space M(s(„) may be identified with the space of nth order differential 
operators with spectral parameter 

i-dtr - Vi{t)i-dtr-^ -... + Vn-2{t)dt - Vn-l{t) - Z. 

For other Lie algebras there also exist canonical representatives for the orbits of the 
gauge action of N+{{t)) on the space of operators of the form (j5.3p (see, e.g., [F3]). 
Recall the element p_i € n_ given by formula (j5.2p . There exists a unique element of 
n+ of the form 

i 

Pi = ^Qei, Ci e C, 
1=1 

such that Pi,p_i, and = form an 5I2 triple. The element € i) then 

defines the principal grading on g such that degpi = l,degp_i = — 1. Let 

be the space of ad -invariants in n, decomposed according to the principal grading. 
Here 

E = {di, ...,di} 

is the set of exponents of g. Then p^ spans V"can,i- Choose a linear generator pj 
of V^an.dj (if the multiplicity of dj is greater than one, which happens only in the 

case g = D^2n,dj = 2n, then we choose linearly independent vectors in Kan.dj)- The 
following result is due to Drinfeld and Sokolov [DSJ (the proof is reproduced in |F2j . 
Lemma 2.1). 

"'^^Note that the action of N+{{t)), and hence the quotient, do not depend on whether we include eez 
in p_i or not. 
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Lemma 4. The gauge action of N^([t)) on the space M(g) is free, and each gauge 
equivalence class contains a unique operator of the form dt + P-i + v(t), where v(t) € 
V^an((^)), so that we can write 



(5.4) v(t) = ^^;,(t).p„ v,{t)(^C{{t)). 

i=i 

Thus, each pomt of the reduced phase space M(3) of the Drinfeld-Sokolov reduction 
is canonicaUy represented by an operator dt + p~i + v(t), where v(t) is of the form 
(ET 



5.2. Classical Hamiltonians. The Poisson commuting Hamiltonians on the reduced 
phase space are constructed in the same way as in the case of regular semi-simple x'- by 
expanding the monodromy of the operators ()5.3p with respect to the spectral parameter 
z. Namely, let M(z) € G be the monodromy matrix of the operator (jS.ip . For any 
invariant function ip on the group G the corresponding function Hip{z) = ip{M{z)) on 
is invariant under the gauge action of N^([t)) and hence gives rise to a well-defined 
function on the quotient M(0). We obtain Poisson commuting Hamiltonians on M(0) 
by expanding H^{z) in z and z^^ . 

The asymptotic expansion at z = oo yields the local Hamiltonians, which generate 
the 0-KdV hierarchy of commuting Hamiltonian flows on the Poisson manifold M(0). 
These Hamiltonians have the form 



Hn = j Pn{vj{t),v'j{t),...)dt, n = di + Nh,N £Z+, 

where di E is an exponent of q and h is the Coxeter number. The integrand P„ is a 
differential polynomial of degree n + 1, where we set degVj^^ = dj + m + 1. 

On the other hand, the z-expansion at z = yields Poisson commuting non-local 
Hamiltonians. It is instructive to consider the z-linear term in the expansion. For that 
we apply Lemma [1] with e = z, 

^^°^w = E/*+^(*)' v{t)€b+{{t)), 

i=l 

and A^^\t) = cq. We obtain that the z-linear term in the expansion of Lp{M{z)), where 
if G (FunG)*^, is given by the formula 

(5.5) i/«= / {ee,^^{t))dt. 



This is the lowest weight component of the single- valued solution ^ip{t) € g of the 
equation 

(5.6) {dt + Ylji + v{t)^^{t) = Q, 

with values in the adjoint representation of g, corresponding to the invariant function 

V9. 
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The commutativity (13. 5p of the Hamiltonians H^p then imphes that 

{<^<} = o 

for all ^ G (FunG)*^. 

Higher order terms in the z-expansion of (p{M{z)) may be expressed as multiple 
integrals of the lowest weight components of solutions of the equation (j5.6p . 

For example, if g = 5(2, the equation (j5.6p in the adjoint representation is equivalent 
to the following scalar third order differential equation: 

(5.7) (5f - 4.vdt - 2v')-^ = 

(its solutions are products of solutions of the second order equation {d^ — v)ij) = 0). It 
has a unique (up to a scalar) single- valued solution ^ip{t) (here we may choose as ip 
the trace function on 5^2, as it generates the algebra of invariant functions on SL2), 
and H^^^ = f ■^^{t)dt. 

Remark 2. A convenient way to compute the higher order terms in the z-expansion of 
^p{M{z)) is to realize the variables of the KdV hierarchy in terms of the variables of 
the modified KdV (mKdV) hierarchy. This provides a kind of "free field realization", 
also known as the Miura transformation, for the commuting Hamiltonians. 
Consider the space MCg) of operators of the form 

(5.8) 9f+p_i + u(t), M{t)eW))- 

The natural map M(g) ^(g) given by the composition of the inclusion M(g) — > M(g) 
and the projection M(g) — > M(g) is the Miura transformation^ It is a Poisson map 
with respect to the Heisenberg-Poisson structure on M(g). 

Because the operator ()5.8p has such a simple structure, it is easier to compute the 
monodromy matrix M{z) for it (and hence for (p{M{z))) than for the operators of the 
form ()5.3p or (j5.4[) . One can then write down more explicit integral formulas for the 
coefficients in the z-expansion in ip{M{z)). For example, in the case when g = sl2 the 
operator ()5.8p has the form 

/ u{t) z \ 
+ V 1 -u{t)) ■ 

The coefficients in the z-expansion of the trace of the monodromy of this operator are 
written down explicitly in [BLZ1| . They are given by multiple integrals of exp(ib2(/)(t)), 
where (j){t) is the anti-derivative of u{t), that is, u{t) = (p'{t). These are the classical 
analogues of the two screening currents discussed at the end of Section 13. 3[ 

Similar formulas may be obtained for other simple Lie algebras. □ 



We have already encountered Miura transformation above, but in a totally different context: 
namely, it appeared in our study of the eigenvalues of the quantum Hamiltonians. Here it appears in 
the study of the phase space on which the classical Hamiltonians are defined. 
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5.3. Quantum Hamiltonians. Now we consider the quantization of the Hamiltonians 
of the g-KdV system. The problem of quantization of the local Hamiltonians has been 
considered, in particular, in [Z | lEYl iKMl IFF41, IFF5| . It is closely related to the problem 
of quantization of the mKdV Hamiltonians and integrals of motion of affine Toda field 
theories (which are in turn related to deformations of conformal field theories). It was 
proved in |FF4l IFF5] that all classical local g-KdV Hamiltonians may be quantized. 
The quantum local Hamiltonians generate a commutative subalgebra in the quantum 
W-algebra ^^(g), where (3 is the parameter related to the central charge. This is an 
associative algebra quantizing the classical (Poisson) W-algebra of local functionals on 
M(g). It is obtained by the quantum Drinfeld-Sokolov reduction of a completion of the 
enveloping algebra of g and depends on one complex parameter, the central charge c, 
which is a function of the level of g (see |FF31 [FB] ). For example, 'W^(s[2) is nothing 
but the completed enveloping algebra of the Virasoro algebra. We will discuss the 
quantum local KdV Hamiltonians in more detail in Section 16.41 below. 

The quantum non-local Hamiltonians were introduced in [BLZl] in the case of g = SI2 
and in [BHK] for g = 5(3 (this construction may be generalized to other g |FFSj ). 

To give an idea as to what these quantum Hamiltonians look like, we discuss the 
quantization of the Hamiltonians H^'^ introduced above. Consider first the case when 

g =3(2. In this case H^p"* is given by the integral of a single- valued solutions ^(^(t) of 
the equation (j5.7p . which is unique up to a scalar. The quantum analogues of solutions 
of this equations are the primary fields <I>(x,3)(?^) of the Virasoro algebra. There are 
actually three primary fields, (^^-^^-^{w) and <^^^^^^{w)^ acting between Verma modules 
over the Virasoro algebra, but only one of them, 3)(^)> ^-^ts from a Verma module 
to itself. This field is the quantization of the unique (up to a scalar) single-valued 
solution ^ip{t) of the equation (j5.7p . Therefore the quantization of the Hamiltonian 

h'^'^ is given by the 0th Fourier component of the primary field <I>^^3^(w). (Here we 
multiply the primary field by ty^t^-^) , where A(x,3) is the conformal dimension, so that 
3)(^) ^ formal power series in integer powers of w.) 

To construct the quantizations of other coefficients in the z-expansion of ip{M{z)), 
we may use the free field realization of the Virasoro algebra in terms of the Heisenberg 
algebra, quantizing the Miura transformation map between the Poisson algebras of 
functions on M(5[2) and M(s[2) described in Remark [2l Then we need to quantize the 
coefficients in the z-expansion of the monodromy matrix of the operator (jS.Sp . The 
resulting quantum non-local Hamiltonians will be given by multiple integrals of the 
screening currents quantizing the classical screening currents exp(ib20(i)) introduced 
in Remark [21 This is the approach taken in |BLZ1| . 

More precisely, the integrals of the two screening currents satisfy the g-Serre relations 
of C/g(s(2), where q = exp(7ri/32) (see (FF5]). The non-local Hamiltonians correspond to 
singular vectors in the Verma modules over C/g(s(2) of critical level (see the discussion 

in Section IS.Sp . In particular, the simplest of these quantum Hamiltonians, H'i^\ is 
given by the g-commutator of the two screening charges. 
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Alternatively, we may use the primary fields ^-^ (w) and g-^ (w), and the fact that 

their integrals generate a quantum deformation of the algebra U{5i2-), as discussed in 
Section 13.31 

This construction generalizes to other Lie algebras. For each g there are primary 
fields <&(i,adj)(^) foi' the corresponding quantum W-algebra 'W^(g), which are the ana- 
logues of the primary fields $(1^3) (w). The vector space of this fields is isomorphic to g, 
and it contains an ^-dimensional subspace (isomorphic to C 3) of those primary fields 
$(1 adj),p(^) which act from each Verma module over 'W^(g) to itself. These are the 
quantum analogues of the single- valued solutions ^'^ (t) of the equation (j5.6p . The 0th 
Fourier coefficients of ^(i^s^di),p{'w),p = 1, ■ ■ ■ ,i (multiplied by tfj^CL^^dj)) g^j-g ^j^g quan- 
tizations of the classical Hamiltonians h!^^ given by formula (j5.5p . Other quantum 
non-local Hamiltonians may be expressed as multiple integrals of the primary fields 
'^(i,adj)(^)- Alternatively, we may use the free field realization of 'W^(g) and express 
then as multiple integrals of the corresponding screening currents [BHKl iFFSj . Note 
that the quantum non-local Hamiltonians commute with each other and with the local 
ones. We will discuss this in more detail in |FF7j . 

The upshot of the above discussion is that we have a large commutative algebra acting 
on the Verma modules (and other highest weight modules) over the Virasoro algebra. 
This algebra is generated by the local and non-local quantum KdV Hamiltonians. We 
now look for a parameterization of the common eigenvalues of these Hamiltonians. 
Our experience with the shift of argument models in which the shift x is regular semi- 
simple suggests that those should be parameterized by affine ^g-opers of special kind 
(discussed in Section [13]) . 

In the case at hand we have a nilpotent shift of argument. In order to understand 
better what kind of opers correspond to this shift, we first consider the analogous 
question in the finite-dimensional case. 

5.4. Nilpotent shift in the finite-dimensional case. Recall from jFFTj that the 
shift of argument Gaudin algebra A^^^ may be obtained as a quotient of the algebra 
Endg 

Ii,x,Kc) where 

\x = I^d^g[[t]]eci ^x- 

Here C,^ is the one-dimensional module over tg[[t]] CI, on which tg[[i]] acts via the 
character 

and the central element 1 act as — /i^. Let us choose x to be the linear functional on g 
taking a non-zero value on eg and zero on all other root generators. The centralizer of 
this X is the Lie subalgebra n_ C g. We wish to understand the corresponding reduced 
algebra (see Section 12. 2p , where ip is a principal character of n_ . 

Hence we perform the quantum Drinfeld-Sokolov reduction of the module Ii^,^ with 
respect to the Lie algebra n_((t)) and its character 



(5.9) fa r 



a = ai,n = 0, 
otherwise 
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(see [FFT]). 

As the result, we obtain a module over the center Z(q) at the critical level, which we 
will denote by //ds(Ii,x)- Recall that Z(g) is a completion of the polynomial algebra 
with generators Si^(^n),i = 1, n G Z, where we use the convention, in which the 

generating currents are 

(the first of then, Si{z) is the unnormalized Segal-Sugawara current), see, e.g., [FFT] . 
It is easy to see that the elements Si^i^m) , i = 1, . . . , — 1; > di, as well as Si^i^m) > iT'e > 
di + 1, act by on the generating vector of Ii^^^ and hence on i?Ds(Ii,x)- Se,(de+i) 
contains a monomial which is the product oi fi l,i = 1,. . . ,i, and eg ®t, with a 
non-zero coefficient. Each factor of this monomial acts by a non-zero scalar on the 
generating vector; for fi®l,i = 1,...,^", this follows from formula (15. Oh . and for 
ee®t this follows from the definition of Ii^^^-. Therefore the monomial itself, and hence 
Si.{di+i)i act by a non-zero scalar on and on //ds(Ii,x)- 

Recall (see, e.g., [F3]) that we have canonical representatives for ^g-opers on the 
punctured disc oi the fornj^ 

where 

t 

(5.10) v(z) =Y,v,{z)-p^, v,{z) G Ciz)). 

i=i 

For c G C, let 0\)"^'^{D) be the space of ^g-opers on the punctured disc of the form 
(ISTOll . where Vj{z) G z^'^^-^<C[[z]] for j = 1, 1, and 

viiz) = cz-'^'-'^ + .... 

Thus, for c = this is the space of ^g-opers with regular singularity, and for c ^ 
this is the space of opers with the mildest possible irregular singularity at the origin. 

Lemma 5. We have an isomorphism 

//Ds(Ii,x)-FunOpi?:;'^(Z5) 

for some non-zero value of c. 

The value c depends on the normalization of p^ as well as the normalization of Xi 
but it is inessential, because by rescaling the coordinate z we may identify the spaces 
OPig'^(-^) foi' all non-zero values of c. 

This has the following immediate corollary: 



Note that in Section [5. II have encountered g-opers as points of the phase space M(0) of the classical 
0-KdV hierarchy. On the other hand, now we are discussing the Hamiltonians in the quantum shift of 
argument Gaudin model associated to g. The spectra of those Hamiltonians are expressed by ^g-opers. 
It is important to realize that these are two different spaces of opers! For more on this, see Section [5.61 
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Corollary 1. For any c € the algebra Ay^ is isomorphic to the algebra of functions 
on the space of ^Q-opers on with regular singularity at G and such that their 
restriction to the punctured disc at oo belongs to Op'^^''^(-Doo)- 

Note that according to the definition of Af^ (see Section 12. 2p , it is isomorphic to 
the center Z{q) of U{q). In addition, the algebra ^lo corresponding to x = is also 
isomorphic to Z{q). Thus, we obtain that SpecZ(0) is isomorphic to Op^l^'^{D) for all 
c G C. ° 

We are now ready to give a conjectural description of the spectra of the KdV Hamil- 
tonians. 



5.5. Spectra of the quantum KdV Hamiltonians. We start with the case of g = 
5l2, which corresponds to the KdV system. 

Recall from Section [4.41 that in the case of regular semi-simple shift x G f) we expect, 
by analogy with the finite-dimensional Gaudin models, that the spectra of the quantum 
Hamiltonians are encoded by the second order differential operators on P^ with global 
coordinate z, which have regular singularity at z = with the leading term —i{£-\-l)/z'^ 
and irregular singularity at oo with the leading term — x^/s^? where s = is a local 
coordinate at oo G P^. These opers are (generically) represented by formula (j4.27p . 

We stress again that this P^ is a curve in the spectral parameter of the affine Gaudin 
model. 

Now, Corollary [1] suggests that for the KdV system, in which x is nilpotent, the 
relevant space of s[2-opers consists of the second order differential operators on P^ that 
have regular singularity at z = with the leading term —i{£ + and irregular 

singularity at oo, but now of order 3, that is, of the form c/s^, c G . By rescaling the 
coordinate z on P^, we can make c = 1, so without loss of generality we will consider 
operators with singularity at oo. The corresponding 512-opers should therefore 
have the form 

where 

(5.11) v{z) 

z~ z \ ' — : / - — : 1^ — Wi)- 

and the coefficients Vj^k in the expansion of v{z) in z — Wj satisfy the equations (j4.19p . 
1 f k \ ^ k 

(5.12) - — Vj^o + Vj^i = 0, j = l,...,m. 



k 1 




4 \w, 



As explained in Section [4.61 this is the condition of "no monodromy" at the points wj. 
As in Section [4.4^ applying gauge transformation by z^'^, we obtain the operator 

(5.13) ^(^)+^^' 
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which we rewrite as the following second order differential operator with spectral pa- 
rameter: 



(5,14) a^-i\r-±±]-^Jl±ll-±-^-f:±^-^. 



^Jt 

W-i I — ^ (z — Wj)"^ ^ Wj z — 



The equations (j5.12p are equivalent to the condition that this operator has no mon- 
odromy around Wj,j = l,...,m, that is, no monodromy on P^, except around the 
points and oo, for all values of A. 

According to our proposal, the differential operators of this kind should encode the 
common eigenvalues of the quantum KdV Hamiltonians on irreducible highest weight 
modules over the Virasoro algebra. Let us recall that the Virasoro algebra obtained by 
the Drinfeld-Sokolov reduction (with respect to the loop algebra of n_) from the affine 
Kac-Moody algebra of level k has the central charge 

(5.15) = 1 - Sttil! 

^ ^ k + 2 

(see |FF1| ). Further, the Drinfeld-Sokolov reduction of the irreducible s[2-module with 
highest weight X = 2i and level k is either (which happens if and only if A; — A € Z+), 
or is the irreducible module over the Virasoro algebra with highest weight (by which 
we mean the eigenvalue of Lq on the highest weight vector) 

A (2£+l)^-(fc + l)^ 

(see [FFT] . p. 316). 

Let us denote by La,c the irreducible module over the Virasoro algebra with highest 
weight A and central charge c (this is the irreducible quotient of the Verma module 
with this highest weight). 

Our proposal is then the following: 

Conjecture 4. The generic common eigenvalues of the quantum KdV Hamiltonians 
on the subspace of Lcf,^At ^ of conformal dimension A^ ^ + m are in bijection with the 
set of operators of the form (I5.14P satisfying the equations (I5.12P . 

For special values of i and k the positions of some of the poles Wj coalesce or approach 
0, and these opers may degenerate, as in the case of regular semi-simple x discussed 
in Section 14. 4i Namely, there may appear points on at which the relative positions 
of the two reductions in the corresponding parabolic Miura oper are given by elements 
of W^\Waff other than sq. If the relative position at a point Wj is given by the element 
so{siSo)"' G VFyVFaff, then the expansion of the function v{z) appearing in (15.14p around 
Wj will have the form 

(2n + l)(2n + 2) , 

^(^) 7 ^2 + • • • • 

[z - wjy 

The coefficients in the expansion of this function at wj will also have to satisfy an 
analogue of equation (j5.12p . which means that our second order operator has no mon- 
odromy at Wj for all A. 
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In addition, when the Verma module Ma^^Ae becomes reducible, some of the points 
Wj may tend to 0. The common eigenvalues encoded by the corresponding limits of 
(I5.14P should then correspond to the eigenvectors which appear in a submodule of this 
Verma module generated by singular vectors. 

It would be interesting to develop an analogue of the Bethe Ansatz procedure for 
constructing eigenvectors of the quantum KdV Hamiltonians. It is not obvious how 
to do this. As we have seen above, the Bethe Ansatz appears when we represent 
our opers as the Miura transformation of Cartan connections. This is possible in the 
Gaudin models associated to regular semi-simple as discussed in Section 14. 4[ But 
the opers ()5.14p cannot be represented as Miura transformation of Cartan connections 
defined globally on P^. (This becomes clear when we look at the restriction of ()5.14p 
to the punctured disc at the point z = oo: the operator has the form + v{s), where 
v{s) = c/s^ + . . ., c E C^, and s = z~^ is the local coordinate at oo; this v{s) cannot 
be represented in the Miura form u(s)^ + dsu{z), where u{s) is a Laurent power series 
in s, as easily checked by substitution.) 

5.6. Two different spaces of opers. In the study of the KdV system, second order 
differential operators, which we call sl2-opers, appear in two different ways. First, they 
appear as points of the phase space MCg) of the classical KdV system, as described in 
Section [5.11 These operators have the form 

df-vit)-z, 

where t is a coordinate on a circle, or a punctured disc, and v{t) is a function on it. 
(Despite what the notation may suggest, t is the "space", not the "time", coordinate of 
the KdV equation.) On the other hand, z is the spectral parameter. The classical KdV 
Hamiltonians are constructed by expanding the monodromy matrix of this operator, 
considered as function of z, near z = and z = oo. 

The s[2-opers appear again as the objects encoding the eigenvalues of the quantum 
KdV Hamiltonians. They have the form 

9^ -t.(z) - Az^ 

where z is the spectral parameter of the classical KdV system. Thus, z is a coordinate 
on P^, and v{z) is a meromorphic function on this P^ (see formula (j5.14p ). On the 
other hand A, is a new spectral parameter. 

It is important to realize that these two spaces of opers are not related to each other. 
In fact, they correspond to two dual theories, one of which is classical and the other 
is quantum. In the classical theory the affine Kac~Moody algebra 5I2 is the central 
extension of 5(2 ((z)), where z is the spectral parameter, and the s[2-opers "live" on a 
circle with coordinate t, which is the space coordinate of the KdV hierarchy. In the 
quantum setting, 5(2 is the central extension of s[2((A)), and the s[2-opers "live" on a 
projective line P^, whose coordinate z is the spectral parameter of the classical system. 

Even more important is the fact that if we replace the Lie algebra 3(2 of the classical 
theory by a more general affine Kac-Moody algebra g, then the corresponding Lie 
algebra of the quantum theory will be not g, but its Langlands dual affine Lie algebra 
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^g, as we will discuss below. This makes it more clear that the classical and quantum 
theories are really different. 

5.7. Comparison with the conjecture of Bazhanov, Lukyanov, and Zamolod- 
chikov. We wish to compare our proposal to the conjectural relation between the 
eigenvalues of the quantum KdV Hamiltonians and second order differential operators 
on of special kind, proposed by Bazhanov, Lukyanov, and Zamolodchikov in |BLZ5j 
(generalizing the proposal made earlier in [DTll IBLZ4j for eigenvalues on the highest 
weight vectors). Understanding the underlying reasons for this relation was in fact one 
of the main motivations for this paper. 
Let us assume that k ^ —2 and set 

k + l 
a = —- . 

k + 2 

We will apply the change of variables 



(2a + 2)2 

to the differential operator ()5.14p . 

The general transformation formula for a second order operator 

52 - v{z) : 17-1/2 ^ ^3/2 
under the change of variables z = ip{x) is 

2 1 

(5.17) v{z) ^ v{ip{x)) {if') --{(p,x}, 

where 

99'" 3/99"' 2 



(5.18) W,x] , , 

is the Schwarzian derivative of ip\^ 

After a straightforward calculation, we find that the new operator looks as follows: 

(5.19) dl - - + 2-^f2 log(^'"'-' - ^j) + ^> 
where 

+ 1) = 4(a + 1)2^(^ + 1) + ^2 + 2a + ^ = 4(a + l)Afc,^ + " ^> 
Zj = (2a + 2fwj, 
E = -(2a + 2)^A 

This operator coincides with the Schrodinger operators parameterizing the spectra 
of the quantum KdV Hamiltonians on i^c^.A^fe introduced in jBLZ5] (formula (1))E1 



-'^^Note that it is imperative to consider our second order operators as acting from Q,^^^'^ to Q?'^'^ , 
for this is the only way to ensure that their property of having the principal symbol 1 and subprincipal 
symbol is coordinate-independent. 

-'^^Note that i corresponds to the i in [BLZ5] , but it is not equal to our £. 
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Moreover, our condition ()5.12p means that the operator (I5.19P has no monodromy 
around the points Zj for all E and is therefore equivalent to the algebraic equations of 
[BLZ5] (formula (3)). 

Thus, we see that our proposal reproduces the proposal of [BLZ5] in the case of the 
quantum KdV system! We view this as a confirmation of the correspondence between 
the spectra of the Hamiltonians of quantum soliton systems and opers for the Langlands 
dual affine algebras which we propose in this paper. 

5.8. Generalization to the case of 5[„. It is easy to generalize formula (I5.14p to the 
case of sin- In this case the Langlands dual Lie algebra is also sin- The s(n-opers may 
be represented as nth order differential operators with spectral parameter. We wish 
to describe the differential operators of this kind on that encode the spectra of the 
quantum Hamiltonians of the 5[,i-KdV system acting on the irreducible module over 
the W-algebra 'Wfc(s(„) obtained by the quantum Drinfeld-Sokolov reduction of the 
irreducible g-module L^^^ of highest weight € [)* and level k. The quantum Drinfeld- 
Sokolov reduction is described in [FF31 IFKWj . In particular, the central charge of 
Wfc(s[„) is given by the standard formula 

c=l + 2n{n+l)-n{n+l)[{k + n)+ ^ 



k + 

These operators have the following general form: 

(5.20) [-d,T - Vl{z){-d,T-^ - . . . + Vn-2{z)d, - Vn-l{z) - Z^A, 

and act from fi^^""^^/^ to fi^"'"''-'^)/^. This determines their transformation properties 
under the changes of coordinate z; in particular, this ensures that the property that 
their principal symbol is (—1)" and the subprincipal symbol is is preserved by the 
changes of coordinate. The coefficients Vr{z) will be rational functions on with poles 
at z = 0, (X), and finitely many other points Wj,j = l,...,m (where m is the degree of 
the corresponding module over Wk^sin))- These coefficients should satisfy the following 
properties: 

(1) At 2 = oo the operator (I5.20p has the mildest possible irregular singularity; namely, 
we have 

Vr{s) '-^ + - - - , r = l,...,n-2; 



s 



Vn-lis) ^ + - - - , 

where Vr{s) are the coefficients of the operator obtained from (I5.20|) by the change of 
variables z s = z~^. Here, as in the case of si2, we normalize the leading coefficient 
of Vn-i{s) to be equal to 1. This may always be achieved by rescaling z. 
(2) Near z = the operator (I5.20p has regular singularity, that is 

where the coefficients Cr(j^) are determined by the highest weight i' of the g-module 
Lj^^fc. Namely, representing u as (z^i, . . . , z/„), where G C and Y17=i ~ we find 
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the Cr{i^)'s from the fohowing formula: 



(5.21) i-d.r - E "-^(-d^r-'-' = (-d. + ^) . . . (-a. + 



i=l 



Z 



(3) At the points wj the operator (j5.20|) has regular singularity, 

Cr{9) 

where 9 = (1,0... ,0,-1) is the maximal root of s[„, which is the highest weight of 
the adjoint representation. There is also the additional requirement that the operator 
(j5.20p has trivial monodromy around the point Wj for each j = 1, . . . ,m, and all A. 

Now, our proposal is that the nth order differential operators of this kind should 
correspond to the generic common eigenvalues of the quantum s[„-KdV Hamiltonians 
on the degree A + m subspace of the irreducible module obtained by the quantum 
Drinfeld-Sokolov reduction of the irreducible g-module Lj^ ^ of highest weight u £ i)* 
and level k. In general, some of the poles wj may coalesce, as in the case of 5(2. 

On the other hand, in |BHK| third order differential operators have been associated 
to the common eigenvalues of the quantum sta-KdV (Boussinesq) Hamiltonians on the 
highest weight vectors in irreducible representations of the W-algebra associated to sts. 
Analogous formulas have been proposed in [DDTlt IDMSTl IDDT2j in the case of s[„. 
As in the case of 5l2 discussed in the previous section, we expect that by applying a 
change of variables to the above operators when the set of points Wj,j = 1, . . . ,m, is 
empty (this corresponds to the highest weight vector), we will obtain the operators of 
[BHKi IDDTll IDMSTl IDDT2j . For m > we obtain a generalization of the proposal 
of [BHKI IDDTll IDMSTl IDDT2j to non-highest weight vectors. 

Explicitly, this change of variables is given by the formula 

Z = — 

(na + n)" 

where 

k + n — 1 
a = ; . 

k + n 

It is easy to see that under this change of variables the irregular term of our differential 
operator gives rise to the term x"'°', and the term Ax'^ gives rise to the new spectral 
parameter term 

no. 

E = {na + n) "+i A, 

which is independent of x. The poles of the new operator will be at the points x = 0, 00 
and Zj,j = l,...,m, where 

Zj = [na + n)^Wj. 

If the affine Lie algebra is of classical type, we may also realize affine opers as 
scalar (pseudo) differential operators, following [DSJ. It is possible to describe those of 
them that encode the spectra of the quantum Hamiltonians of the g-KdV system in 
a similar way to the case of sin- We expect that the corresponding operators (in the 
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special case when the set of extra points wj is empty) are related to those proposed 
in [DDTll IDMSTl IDDT2| after a change of variables. However, if we wish to describe 
these affine opers for an arbitrary ^g, we should use instead the canonical form for 
opers discussed in the next section. 

5.9. Generalization to an arbitrary afRne algebra g. It is also easy to generalize 
formula (15.14P to the case of an arbitrary affine Kac-Moody algebra. For simplicity 
we will assume here that is an untwisted affine Kac-Moody algebra I, where I is a 
simple Lie algebra, but the general case may be analyzed along the same lines. 
Let us recall (see, e.g., [F2j ) that [-opers have canonical form 

e 

V = + p_i + ^ Vriz) ■ Pr, 
r=l 

where Pj.,r = are the elements defined in Section [5.11 Here vi{z) transforms 

as a projective connection (that is, by formula (j5.17p ) and Vr{z),r > 1, transforms as 
a {dr + l)-differential under the changes of coordinate z. Our l-opers on will have 
the form 

e 

(5.22) d,+p_^ + Y^Vr{z)-p^ + Xz''pi, 

r=l 

where each Vr{z) is a rational function on with poles at z = 0, oo, and finitely many 
other points Wj,j = 1, . . . , m. These coefficients should satisfy the following properties: 

(1) At z = oo the operator ()5.22p has the form 

^-(*)~^ + ---' r = !,...,£- 1; 

^ , , 1 

where Vr{s) are the coefficients of the operator obtained from (|5.22|) by the change of 
variables z s = z~^ . 

(2) Near z = the operator ()5.22p has regular singularity, that is 

where the coefficients Crii') are determined by the highest weight v of the g-module 
-Ljy_fc. They are determined by the following rule: the element 



P-l + (criu) + ^5r,l] Pr 

r=l ^ ^ 



is the unique element in the Kostant slice of regular elements, 

p_i+v, V G b, 
which is conjugate to p_]^ — z/ (see [F6j . Sect. 9.1). 
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(3) At the points Wj the operator ()5.22|) has regular singularity, 

where 6 is the maximal coroot of I. In addition, we require that the oper have trivial 
monodromy around the point Wj for each j = 1, . . . , m, and all A. 

Our proposal is that the ^g-opers (|5.22|) satisfying these conditions encode the generic 
common eigenvalues of the quantum Hamiltonians of the g-KdV system on the irre- 
ducible modules of the corresponding W-algebra obtained by the quantum Drinfeld- 
Sokolov reduction \FF3\ IFKWj from the irreducible module over g with highest weight 
ly and level k. 

6. Affine Gaudin models with regular singularities 

In the previous section we have discussed the affine Gaudin models corresponding 
to the shift of argument method. We have seen that the classical Hamiltonians of 
these models are integrals of motion of familiar soliton hierarchies, such as the KdV or 
AKNS hierarchies. Hence the quantization of these Gaudin models naturally leads to 
the corresponding quantum Hamiltonians and their diagonalization. However, in the 
universe of all Gaudin models (see |FFT| ). the shift of argument Gaudin models are 
very special: they correspond to two points on P^, one with regular singularity and one 
with irregular singularity of order 2. It is natural to develop a similar theory for more 
general, multi-point, Gaudin models. 

The simplest of these models are the ones with "regular singularities" at finitely 
many points of P^, in the terminology of [FFTj . In this section we study the corre- 
sponding affine Gaudin models. We start out by briefly summarizing the main features 
of the Gaudin models with regular singularities in the finite-dimensional case. We then 
discuss analogous patterns in the affine Gaudin models and their potential connection 
to the quantum KdV system. In particular, we will suggest an explicit construction 
of eigenvectors of the quantum KdV Hamiltonians on representations of the unitary 
minimal models via Bethe Ansatz. 

6.1. Gaudin model, finite-dimensional case. Let g be a simple finite-dimensional 
Lie algebra over C. Recall that it carries the invariant inner product kq normalized in 
such a way that the squared length of the maximal root is equal to 2. Let {J"} be 
a basis of g and {Ja} the dual basis with respect to this inner product. The Gaudin 
Hamiltonians are the elements of [/(g)®^ given by the following formula: 

(6.1) -^,^Y.T.-Th-' i = i,...,iv. 

Here = 1 ® . . . ® A® . . . <^ 1, with A in the ith position. 

As shown in [FFR l [F3] . the quadratic Gaudin Hamiltonians are included in a large 
commutative subalgebra of C/(g)®^, called the Gaudin algebra. According to [F3j , 
this algebra is isomorphic to the algebra Fun 0pi^(P-^)(2.) oo of functions on the space 
OPlq (IF'"'^)(2i),oo of ^g-opers on P^ with regular singularities at the points zi.,...,zn 
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and oo. Here is the Langlands dual Lie algebra of g, whose Cartan matrix is the 
transpose of the Cartan matrix of q. 

This isomorphism comes about [FFR^ IF3| in the same way as in the case of the shift 
of argument Gaudin models, discussed in Section [2l Namely, using the coinvariants (or 
conformal blocks) of modules over the affine algebra g of critical level, we identify the 
Gaudin algebra with a quotient of the algebra EndjVcrit, where Vcrit is the vacuum 
module of critical level (see formula (j2.5p ). The latter algebra is isomorphic to the 
algebra of functions on the space OpLg(L>) of ^g-opers on the (formal) disc [FF3l IF4] . 
This implies that the Gaudin algebra is isomorphic to the algebra Fun Opi^(P^)(2.)^oo 
of functions on Op?^^(P^)(^^)^oo- 

As before, the appearance of the Langlands dual Lie algebra ^g here is important. 
It signifies the fact that the above isomorphism may be viewed as a special case of the 
geometric Langlands duality (for more on this, see [Fl^ \F5\ IFFT] ). 

Thus, we obtain that the common eigenvalues of this Gaudin algebra, including 
those of the Gaudin Hamiltonians Hj, are encoded by ^g-opers in 0pi^(P^)(2.) oo- Fo'^ 
example, in the case of g = 5(2 the corresponding opers are the second order operators 
of the form 

^ A. 



(6-2) s'-Er^^-E^^- 

^(.-z.)2 

In particular, the eigenvalues of the Gaudin algebra on the tensor product of the Verma 
modules (^^1 M,y- are encoded by such differential operators, in which 

\ = ^Vi{ui + 2); 

this is the value of the Casimir operator in [/(sb) on the ith factor Mj^-. The operator 
()6.2p is then determined by the remaining parameters Q, i = 1, . . . , A^, and in fact they 
are nothing but the eigenvalues of Hj, i = 1, . . . , (see (F3j ). 

Thus, in the case of g = ^[2 the Gaudin algebra is generated by the Hj's and the 
Casimir operators along each factor in U{5\2)®^ ■ But for all other simple Lie algebras 
there are higher order Gaudin Hamiltonians, whose eigenvalues are recorded by poly- 
nomials in the coefficients appearing in the expansions of the ^g-opers near the points 
Zi. For instance, for g = s(„ the role of (16. 2p is played by the nth order differential 
operators with regular singularities: 



where Vj{z) has the form 

N j+1 



= EE' 



C-ijk 



.=1 k=l - 

Thus, such an operator is determined by the coefficients Cijk appearing in its expan- 
sion around the singular points = 1, . . . , A. The quadratic Gaudin Hamiltonians 
correspond to Q^i^i, and the coefficients with j > 1 correspond to higher Gaudin Hamil- 
tonians of order j + 
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6.2. Bethe Ansatz. Let us choose a cohection of weights vi,. . . ,V]<[ S f)*. Then we 
have the Verma modules My^ with highest weights = 1, . . . ,un- Eigenvectors of 
the higher Gaudin Hamiltonians on the tensor product (^j^^ My^ may be constructed 
by Bethe Ansatz (see [FFRl [F3] ). 

Namely, for any collection of distinct complex numbers wi, . . . , Wm such that wj ^ Zi 
for all i and j, and a collection ctj^ , • • • , on^ of simple roots of Q (so that ij € I, the set 
of vertices of the Dynkin diagram of g) consider the following Bethe vector 



(6.3) <^(<,...,u;j;r 



, TTt \r — ^ r^^T (g)...u^^. 

, , „s , , (^ife ~ W,k)[W,k - W,k) . . . iw,k - Zk) 

p=(pi,...,P^) fc=l ^1 •'a I'i' ^ J^k 

(k) 

Here denotes the generator of the Lie algebra n_ C g corresponding to the ith 
simple root, acting on the kth. factor of ^jLi ^uj^ ^-iid the summation is taken over all 
ordered partitions P-*^ U U . . . U of the set { 1 , . . . , m} , where = { , , . . . , ja^. } 
and {ii, «2) • • • ) ^aj.} ^'^^ the corresponding labels in /. 

It is proved in |FFR1 IF3j that the vector (^{w^i , ■ ■ ■ , w^m) is a common eigenvector of 
all higher Gaudin Hamiltonians if the following system of equations is satisfied 

^ Wj - Zi ^ Wj - Ws 

where (•, •) denotes the pairing between f)* and f). These are the Bethe Ansatz equations 
of the Gaudin model. 

As shown in [FFRl IF3j , the common eigenvalues of the higher Gaudin Hamiltonians 
on the Bethe vector 4>{wY , . . . ) given by formula ()6.3p (provided that it is non- 
zero) are encoded by the ^g-oper obtained by applying the Miura transformation to 
the Cartan connection on the ^//-bundle O'' over P^, 



N m 
^-^ Z — Zi ^-^ 

i=l j=l 



Wi 



with regular singularities at the points Zi and Wj. The Miura transformation maps 
such a connection to an ^g-oper on with regular singularities at these points. This 
oper has no singularities at the points Wj if and only if the Bethe Ansatz equations 
(j6.4p are satisfied (see [FFRl IF3j ) . Equivalently, we can phrase this as saying that the 
Bethe Ansatz equations mean that the Miura transformation of the connection V have 
no monodromy around Wj^j = 1, . . . , m. 

There is an additional restriction on the common eigenvalues of the Gaudin Hamil- 
tonians that may occur in the tensor product of finite-dimensional representations, and 
that is the condition of no monodromy around the points Zi,i = 1, . . . ,N, and oo. De- 
note by Vu the irreducible finite-dimensional g-module with highest weight z^, which is 
dominant and integral. The higher Gaudin Hamiltonians act on V^. and commute 

with the diagonal g-action. Let us describe the joint eigenvalues of these Hamiltonians 
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on the isotypical component 




where Voo is another dominant integral weight. According to |F2| . Corohary 4.8 (see 
also |F3J, Theorem 2.7,(3)), the set of ^g-opers representing them (counted without 
multiplicities) is a subset of the set of ^g-opers on with regular singularities at 
Zi,i = 1, . . . ,N, and oo, with fixed residues determined hy Viji = 1, N , and i^ooi and 
with trivial monodromy representation (such opers may equivalently be described as 
z/j-regular at Zi and z^oo-regular at oo, see below). Furthermore, according to Conjecture 
1 of [F3j . this subset is in fact the entire set of such opers. In other words, all opers of 
this kind correspond to common eigenvalues of the Gaudin Hamiltonians. 



6.3. Generalization to the afRne case. Now we would like to generalize the above 
results on the Gaudin models to the case when a finite-dimensional simple Lie algebra 
g is replaced by an affine Kac-Moody algebra g, by analogy to the shift of argument 
models discussed above. We begin by summarizing the known facts about the affine 
Gaudiri^model. 

Let A be the set of roots of an affine Kac-Moody algebra g (see Section [3.ip . The 
Lie algebra g has a canonical invariant inner product kq whose restriction to g (8) 1 is 
equal to kq. Let us choose dual bases {J"} and {Ja} of g with respect to kq. Then we 
choose dual bases {J°^}, {Ja} of g with respect to kq in the following way. The former 
consists of J" t", n € Z, 1 and d, and the latter consists of Ja <^ t^,n (z d and 1. 

Let zi,...,zjsr be a set of distinct complex numbers. The quadratic Gaudin 
Hamiltonians associated to g are given by the same formula (j6.ip as in the finite- 
dimensional case: 

(in fact, such Hamiltonians may be defined for any symmetrizable Kac-Moody algebra). 
The elements belong to a completion of the A^-fold tensor power of the universal 
enveloping algebra f/(g). It gives rise to a well-defined linear operator on the tensor 
product (^j^j^Mj, if each Mi is a smooth g-module. By definition, a g- module M is 
called smooth if for any f G M we have g ® t™'C[[t]]u = for sufficiently large m € 
In particular, highest weight modules, or modules from the category 0, are smooth. 
It is clear that the Hamiltonians Hj commute with the diagonal action of g and that 

N 

(6.5) 5]H, = 0. 

i=l 

For a weight € f}*, denote by Mjy the Verma module 

Mi. = Ind! 
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Consider the problem of diagonalization of the quadratic Gaudin Hamiltonians on the 
tensor product of Verma modules (^^^ M^,^. The Bethe Ansatz presented in the pre- 
vious section generalizes in a straightforward way to the case of afhne algebras (and, 
moreover, to all symmetrizable Kac-Moody algebras). 

More precisely, for any collection of distinct complex numbers wi, . . . , Wm such that 
Wj 7^ Zi for all i and j, and a collection , • • • , ctim of simple roots of g (so that ij G /, 
the set of vertices of the Dynkin diagram of 3) we have the Bethe vector given by formula 
(|6.3p . Then, according to the results of [FFRj . Sect. 5 (based on the computations 
of |SV1| ). this vector is an eigenvector of the quadratic Gaudin Hamiltonians = 
1, . . . , A^, if the Bethe Ansatz equations (j6.4p are satisfied0 This is in complete analogy 
with the finite-dimensional case. 

Thus, so far in the affine case we have the quadratic Gaudin Hamiltonians and 
the Bethe vectors ()6.4p . However, by analogy with the finite-dimensional case (and 
the affine shift of argument Gaudin models), we expect they are included in a large 
commutative subalgebra of a completion of U{q)®^ . This raises the following questions: 

(1) Are there higher order affine Gaudin Hamiltonians? 

(2) What parametrizes the joint eigenvalues of the affine Gaudin Hamiltonians? 

In the next section we will take up the first question in the simplest case of two points 
-21,^2 (or, more properly, three points, since we should count 00 as an extra point). In 
this case there is only one quadratic Gaudin Hamiltonian, Hi, because H2 = —Hi due 
to the identity (16. 5p . We will show that Hi is nothing but the operator Lq of the coset 
Virasoro algebra acting on the tensor product of two g-modules Mi M2. Therefore 
all quantum KdV Hamiltonians H2i+i,i > 0, written with respect to this Virasoro 
algebra, commute with Hi (which is equal to Lq = Hi). We will conjecture they these 
are precisely the higher (local) affine Gaudin Hamiltonians. 

In the case of multiple points we do not have an analogous formula for the lo- 
cal Gaudin Hamiltonians, but we will give a procedure for constructing the non-local 
Gaudin Hamiltonians (by analogy with the above construction of the non-local quan- 
tum AKNS and KdV Hamiltonians). We will also consider the quasi-classical limit of 
this construction, which is an affine analogue of the Hitchin system. 

We will then discuss the second question of describing the spectra of the affine 
Gaudin Hamiltonians. By analogy with the finite-dimensional case, and the affine shift 
of argument Gaudin models, we will conjecture that the joint eigenvalues of the affine 
Gaudin Hamiltonians are parameterized by the affine ^^-opers with regular singularities 
at the points zi, . . . , zjy and 00. 

6.4. The case of two points. Let us consider the simplest case of two distinct points 
zi,Z2. Thus, we have two smooth g-modules Mi and M2 and the quadratic Gaudin 
Hamiltonian 

H = ^Jh0 

a 



'This may also be proved along the lines of [RV] , 
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acting on Mi (g> M2. Note that Si = — H2 = S/(zi — 22)- It commutes with the 
diagonal action of g on Mi M2 . Our question is to construct the higher affine Gaudin 
Hamiltonians, which commute with the diagonal action of 5 and with H. 

We will construct them using the coset construction of [GKOj . Suppose that the lev- 
els ki and k2 of the modules Mi and M2 are non-critical, i.e., ki 7^ — /i^ and in addition 
ki +k2 7^ —hy. Then the GKO construction gives us the Virasoro algebra which com- 
mutes with the diagonal action of g in the completed tensor product t^fc^(g)(8)C/fc2('0), 
where ?7fc(0) = C/(0)/(l — k). This algebra is generated by -L^^^,n € Z, such that 

nGZ 

Here T^^(w) is the Segal-Sugawara current at level k 7^ — /i^, 

and the lower indices refer to the first, the second or the diagonal g in C/fci (g)®f/fc2(s)- 
Without loss of generality we may, and will, assume that the action of Lq^ on the 
g- modules Mi and M2 coincides with the action of — d G g. Then we find that 

-{ki + k2 + /i^)L[f^° = E E -^"'^ ®Jn-ki® Lf - Lf ® k2 

'jz0r = E. 



Thus, up to a non-zero scalar, the quadratic Gaudin Hamiltonian coincides with the 

generator L^^'~' of the coset Virasoro algebra. Since the entire coset Virasoro algebra 

commutes with the diagonal action of g, the first requirement on the higher Gaudin 

Hamiltonians is satisfied for any element of its completed universal enveloping algebra. 

Now we search for a commutative subalgebra of this enveloping algebra containing 
7- GKO 

° ■ 

Consider the case of g = SI2. Then a natural candidate for this commutative sub- 
algebra is provided by the algebra of quantum KdV Hamiltonians. Let us recall that 
the quantum KdV Hamiltonians are certain elements H2i+i,i > 0, of the completed 
universal enveloping algebra of the Virasoro algebra whose symbols are the classical 
KdV Hamiltonians (see Section [5. 2p . The first of them. Hi, is nothing but Lq, and the 
next one is 

^3 = ^^0 + E = 1 [ ■'^i^f- w^dw. 

n>0 '' 

The other Hamiltonians have the form 

H2i+i = j :T{wY^^:uP'^^^dw + \oweT order terms. 

Explicit formulas for the first few of them may be found in [BLZlj . The existence of 
these Hamiltonians has been established in [FF41 [FF5] . We conjecture that 
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The quantum KdV Hamiltonians of the coset Virasoro algebra are the higher 
Gaudin Hamiltonians of the two-point 5I2 affine Gaudin model with N = 2. 

More concretely, we conjecture the following: 

Conjecture 5. The affine Bethe vectors ()6.3p (for sl2 and N = 2) are simultaneous 
eigenvectors of all quantum KdV Hamiltonians of the coset Virasoro algebra, provided 
that the Bethe Ansatz equations (16. 4p are satisfied. 

Consider the special case when the module Mi is the vacuum module V(o,i) of level 
ki = 1 and M2 is an integrable module V(^n,k) of level k2 = k £ and highest weight 

n £ Z+,n < k, with respect to the one-dimensional Cartan subalgebra f) C f). Then, 
according to |GKQ| . we have the decomposition 

fc+i 

^(0,1) ^(n,fe) - ^ ^(m,fc+l) "Xi Ln,m, 

with respect to the diagonal action of sl2, where Ln^m is the irreducible module over 
the Virasoro algebra from the (fc + 2, A; + 3) unitary minimal model with the central 
charge 1 - 6/(/c + 2)(A; + 3). 

If the above conjecture is true, then we obtain an effective method for construct- 
ing eigenvectors of the quantum KdV Hamiltonians on the irreducible modules of the 
Virasoro algebra from the unitary minimal models. 

In the case of a general simply-laced Lie algebra g, we also expect that the higher 
affine Gaudin Hamiltonians coincide with the quantum g-KdV Hamiltonians, which 
generate a commutative subalgebra in the coset W-algebra. If this is true, then the 
Gaudin Bethe vectors would give us eigenvectors of the quantum g-KdV Hamiltonians. 

In the above discussion we considered the "local" Gaudin Hamiltonians, that is, those 
that may be expressed as infinite sums of monomials in the generators of g of bounded 
order (for example, the Hamiltonian H has order 2). However, in addition to these local 
Hamiltonians there are also non-local ones. In order to explain the appearance of these 
non-local Hamiltonians, we consider the quasi-classical analogues of the affine Gaudin 
models: the affine Hitchin systems. 

6.5. AfRne Hitchin systems. The Gaudin models associated to simple Lie algebras 
have quasi-classical versions, which may be interpreted as the Beauville-Hitchin systems 
on with ramification (see [HI IBea^ IMarl IBotj ) . In the case of Gaudin models with 
regular singularities the algebra of quantum Hamiltonians is a commutative subalgebra 
of C/(g)®^. Its quasi-classical counterpart is then a Poisson commutative subalgebra in 
S{q)^^ = Fun((g*)®^). Let us identify the Poisson manifold (g*)®^ (with its Kirillov- 
Kostant Poisson structure) with the space of one-forms on with poles of order one 
at the points 21, ... , zn, and 00, 

N . 

(6.6) L = V -^dz, A, e g*. 

Let Pi, . . . , Pi be generators of the algebra Invg* of G- invariant functions on g*, of 
degrees di + 1, . . . , di + 1, where the dj's are the exponents of g. Then we have the 
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Hitchin map 

TV 

p : ^ :k = ifO {¥\ {n{z, + ... + zj, + oo))«('^»+^) 

i=l 

L^{Pi{L),...,Pe{L)). 

The pull-backs of polynomial functions on Ji to (g*)®^ via the map p give rise to 
Poisson commuting Hamiltonians on (g*)®^. They generate a Poisson commutative 
subalgebra of 5(0)*^^. According to the results of \F3\ IFFTj . this is precisely the 
algebra of symbols of the higher Gaudin Hamiltonians. 

Now we wish to generalize this construction to the affine case, where we will have 
commuting quantum Hamiltonians in a completion of U{q)^^ . Their symbols should 
therefore belong to a completion of the algebra of functions on (g*)®''^. It is natural 
to identify the latter with the space of one- forms L as in (j6.6p . but with each Ai now 
being an element of q*. Actually, as in the case of the shift of argument subalgebra, 
we will consider instead of the ith copy of g*, the hyperplane g^. in g* which consists 
of the linear functionals on g taking value ki on the central element 1. We will identify 
this hyperplane with the space of fcj-connections on the punctured disc 

kidt + A{t), A{t)eg{{t))^cSg*{{t)), 

so that the coadjoint action of G([t)) becomes the gauge action on the fcj-connections 
(compare with Section [32|)- This hyperplane carries a canonical Kirillov-Kostant Pois- 
son structure. 

Then the L-operator of our affine Hitchin system may be represented in the form 

(6.7) L = J2^^^^±^dz, AGQiit))- 

1=1 

In what follows we will assume that ki ^ 0. 

Equivalently, we may rewrite this as a connection on the punctured disc, which 
depends on the spectral parameter z € P^\{zi, . . . , Z]\f, oo}: 

(6.8) L' = dt + j^Y.^dz, 
where 

(6.9) f{z) = ^ 



z 
1 



Now we will define commuting Hamiltonians of the affine Hitchin system by evalu- 
ating on L' gauge invariant functions on the space of connections dt + A{t). 

6.6. Commuting Hamiltonians from the monodromy matrix. As in the case of 
the shift of argument model (see Section 13. 2p , the invariant functions are obtained by 
taking the invariants of the monodromy of the connection with respect to the adjoint 
action. Denote by M{z) S G the formal monodromy matrix of the connection operator 
L' given by formula (16. Sp . As before, M{z) is defined as follows: its action on any 
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finite-dimensional representation y of G is obtained by computing the monodromy of 
the corresponding matrix differential operator dt + A{t) acting on V . 

For any function on G which is invariant under the adjoint action of G, the 
Hamiltonian 

H^{z) = if{M{z)) 

is a gauge invariant function of L' (or of L given by formula (j6.7p ). depending on the 
spectral parameter z. One shows in the same way as in the finite-dimensional case that 

(6.10) {H^{z),H^{y)} = {), \/^,^(^¥nn{Gf, z,yeC\{zi,...,ZN]. 

Therefore we may obtain families of commuting Hamiltonians of the affine Hitchin 
model by expanding the Hamiltonians H^{z) around points of P^. 

Let us first consider the expansion around one of the points Zi,i = 1,. . . ,N. The 
expansion of the operator kiL' at Zi reads as follows: 

k.8. + Mt) (-A.it) + 2 ^ . . , 

where the dots stand for higher order terms in (z — Zj). 

We now apply Lemma[I]in this case, using [z — Zi) as the parameter e. Let us write 

a 

and denote by 

a 

the single- valued solutions of the equation 

(6.11) kidt^i{t) + [Ai{t),-^i{t)]=^, 

corresponding to (/? € (FunG)*^, as before. Then we find, using Lemma [H that the 
(z — Zj)-linear term in the expansion of ip{M{z)), where M{z) is the monodromy of the 
operator fcjL', is equal to 

(0.12) . / E ^ . E ^\ n.m 

(here and below the summation over the repeated index a is understood). The com- 
mutativity ()6.10p of the Hamiltonians then implies that 

{i?i,<^,i/fc,^} = 0, i = l,...,iV, G (FunG)^. 

Higher terms in the expansion of (/9(M(2;)) in (z — Zi) are given by more complicated 
formulas, involving multiple integrals of solutions of (|6.1ip . All of these Hamiltonians 
are non-local, in the sense that they cannot be realized as integrals of differential 
polynomials in the Ai(t). 

In order to obtain local commuting Hamiltonians, we need to consider the expansion 
of the operator V at the points z = yi, . . . , which satisfy /(?/i) = 0, where 

/(z) is given by formula (|6.9p (this is similar to what happens in the more familiar 
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integrable systems; see, e.g., |DSt IRS2j )). We expect that for each of these points 
there is an infinite series of Hamiltonians whose degrees are ah positive integers equal 
to the exponents of g modulo the Coxeter number. These Hamiltonians should define 
interesting hierarchies of commuting flows on our phase space Hi^i 0fc • 

We have not seen these particular hierarchies discussed in the literature, but some 
closely related hierarchies have been studied before. For example, consider the special 
case when N = 2. Then the limit of the L-operator ()6.8p when /c2 — »■ is closely related 
to the L-operator of the principal chiral field model. We expect that the corresponding 
limits of our Hamiltonians coincide with the Hamiltonians of that model constructed 
in [EHMMj (in particular, they have the right degrees). We plan to discuss these 
Hamiltonians in more detail elsewhere. 

6.7. Quantum non-local Hamiltonians. The classical Hamiltonians of the affine 
Hitchin system may be quantized, giving rise to quantum (non-local) Hamiltonians of 
the affine Gaudin model. 

Let us consider first the case of two points, so that N = 2. In this case we have the 
following classical Hamiltonians: 

Hi,^ = j (-^Ji,a(t)+J2,a(t)) n,^mt, 

H2,^ = J (-^J2,a(0 + ^l,a(t)) nA^)dt 

(we have multiplied them by {zi — Z2)). To quantize these Hamiltonians, we replace 
Ji,a{t) by the current Ja\z) and ^i^{t) by the component <^p^'^\z) of the primary field 
<^p\z), introduced in Section [331 

a 

Here, as before, the upper index (?) indicates that this operator acts along the ith factor 
of the tensor product Mi ® M2. We also need to adjust the levels, as in the equation 
(j3.14p (and the Segal-Sugawara formula). 
The resulting quantum Hamiltonians read 

They act on the tensor product M\^^ki <8 of Verma modules over g (or more 

general modules Mi ® M2, provided that the primary fields ^p\w),i = 1,2, are well- 
defined on them). 

It is easy to see that they commute with the diagonal action of g. Indeed, consider 
again the extended affine algebra g = Cd >< g. We choose dual bases { J"^} and { J5} in 
g as in Section lOl Recall that they consist of elements of the form J"'" (resp., Ja,n)-, 
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1 and d (resp., d and 1). Suppose that we have two homomorphisms of g-modules 
Pi : — > Endc Mi, i = 1,2. Then the endomorphism 

(6.13) J^/oi(Ja)®/92(J^) 

a 

of Ml (8) M2 commutes with the diagonal action of q. 

In our case, we choose as the first homomorphism pi the action of q on Mi. Under 
this action 1 maps to ki Id and d maps to the Sugawara operator — Lq- We define the 
second homomorphism as follows: 

1 0. 

One checks, in the same way as for the Segal-Sugawara formula, that this indeed defines 
a homomorphism of g-modules. The resulting operator (j6.13p is the Hamiltonian H2^p- 
Hence it commutes with the diagonal action of q. Since the Hamiltonian Hi^p is obtained 
from -ff2,p by switching the modules Mi and M2, we find that it also commutes with g. 

Furthermore, we claim that the Hamiltonians Hi^p commute with each other for all 
i = 1,2 and p = !,...,£. This may be checked by a direct calculation which will be 
presented in [FF7] . 

The above formulas are generalized in a straightforward way to the case of points. 
We now have the following commuting Hamiltonians quantizing the classical Hamilto- 
nians (pT2]) : 

where z = 1, . . . , A^, p = 1, . . . , They act on the tensor product 0^i Mx.^^i of Verma 
modules, or more general highest weight g-modules of levels ki,i = 1, . . . , N (recall that 
ki 7^ — /i^ by our assumption). These are the simplest non-local quantum affine Gaudin 
Hamiltonians. 

Other non-local quantum Gaudin Hamiltonians may be constructed by generalizing 
the procedure of [BLZll IBHK[ IFFSj , as outlined in Section 13. 3[ We will discuss this in 
more detail in jFF7j . 

6.8. Hamiltonians in the limit Z2 00. It is important to observe that the Hamil- 
tonians of the shift of argument Gaudin model may be obtained from the Hamiltonians 
of the Gaudin model with regular singularities and = 2 in a certain limit as 2:2 ^ cxd 
(we also set zi = 0). In the finite-dimensional case this was explained by L. Rybnikov 
in [Rj. The Hamiltonians of the latter model generate a commutative subalgebra in 
U{q) ® U{q). Consider the limit Z2 — > cjo, in which we also degenerate the second 
factor U{q) (corresponding to the point Z2) to S{q). Then we obtain a commutative 
subalgebra in U{g) S{g). Applying the homomorphism S{g) — > C, corresponding to 
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evaluation at x ^ 0*i we then obtain a commutative subalgebra in U{q). It turns out 
to be the quantum shift of argument subalgebra A^, as shown in [R] . 

We may apply the same limiting procedure in the affine case. The above quantum 
Hamiltonians of the Gaudin model with regular singularities and N = 2 generate a 
commutative subalgebra in the (completed) tensor product U{Q)(dU{Q). Taking the 
limit ^2^00 and at the same time degenerating U{q), we obtain a commutative 
subalgebra in U{q)®S{q). We then need to evaluate our Hamiltonians at x ^ 0* along 
the second factor. For example, let us apply this procedure to the Hamiltonians Hi^p 
and H2^p constructed above. Let us set k2 = 0. Then in our limit the generators Ja,n of 
becomes classical, that is, linear functional on g*. Evaluating them on x G f)* C g*, 
we obtain precisely the non-local Hamiltonians 

a 

introduced in Section [3.31 

On the other hand, if we apply this procedure to H2^p, then the primary fields ^p{w) 
will become classical, i.e., will degenerate to ^'^{w), and we will need to evaluate them 
on X- The resulting Hamiltonians will be linear combinations of Jafi that belong to 
f) C 0, where i) is the centralizer of x in 0- Thus, we obtain local (in fact, linear) 
Hamiltonians of the shift of argument affine Gaudin model. 

Similar formulas will appear in general. The non-local Hamiltonians in the Gaudin 
model with regular singularities and N = 2 have, roughly, the following form: 



^Ai(g)Bi or ^Si^Ai, 



where the Ai are non-local expressions which are quantizations of the derivatives of the 
"quantum monodromy matrix", such as <I>p (-n)' local expressions, that 

is, of finite degree in J^. When we degenerate U{q) to S{q) along the second factor, 
either Bi or Ai become classical, and we may then evaluate them on x G 0*- In the first 
case we obtain non-local Hamiltonians of the shift of argument affine Gaudin model 
(such as the g-AKNS model). In the second case, we obtain local Hamiltonians. 

A similar limiting procedure may be constructed for the Hamiltonians of the quantum 
0-KdV systems. This suggests a possible way to construct local quantum Hamiltonians 
of the g-AKNS and g-KdV systems. 



6.9. Spectra of the quantum Hamiltonians and afRne opers. In |F2j to each 
solution of the Bethe Ansatz (j6.4p equations for the affine Lie algebra g the following 
connection on the ^if-bundle QP was attached: 

N 

(6.15) V = a, - V + 

^-^ Z — Zj 

i=i j=i 



^ Z — Wi 



It was shown in |F2j that the Bethe Ansatz equations are equivalent to the condition 
that the ^g-oper V = V + p_i obtained by applying the Miura transformation to 
the connection (j6.15p is regular at the points Wj,j = l,...,m, or, equivalently, has 
no monodromy around these points. This oper has regular singularities at the points 
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Zi,i = 1,...,N and oo with the residues determined by the highest weights = 



This suggests that the true parameters of the Bethe vectors are not the Cartan 
connections (I6.15p . but the corresponding ^g-opers obtained by applying the Miura 
transformation. Thus, we propose the following answer to question (2) of the previous 
section: what parameterizes common eigenvalues of the affine Gaudin Hamiltonians? 

The common eigenvalues of the affine Gaudin Hamiltonians on the tensor 
product ^iLi M^. of Verma modules over g are encoded by ^g-opers on 
with regular singularities at the points zi, . . . ,zn and oo. 

Again, we emphasize that this is a curve in the spectral parameter z of the affine 
Gaudin model. 

To see what these ^g-opers look like, consider the case when = ^0 = 5(2. We will 
use the conventions of Section 14. 4[ We again divide the set of points wi, . . . , Wm into 
two subsets wl, . . . ,wj^_^, and w^,. . . ,w'^^, corresponding to the simple roots ai and 
ao, respectively. The connection (|6.15p reads 



1,...,N. 




where 




N 



(6.16) 



Applying the gauge transformation by 




1 -u{z) 
1 




we obtain the operator 



(6.17) 




where 



(6.19) 



(6.18) 
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These formulas follow from the Bethe Ansatz equations (16. 4p . which read 

TV ^ 1 ™° 1 

(6.20) T.— T + E^ 0=0' J = !,•••, mi, 

1=1 J s^j * s=\ ] * 

N k,_n, rn-, 

(6.21) Y^^^ + Y^ T-Et 0=0' J = l,---,rno. 

1=1 J s=l J * s^j J * 

These equations imply, in particular, that v{z) has no singularities at the points wj, J = 
1, . . . ,mi. 

In addition, because v{z) is obtained as the Miura transformation of u{z) given by 
formula (j6.16p . the coefficients of its expansion at each point w^j,j = 1, . . . ,mo. 



0^n 



satisfy the equation (j4.19p . This is the condition that the operator (j6.17p has no 
monodromy around the point w^y For instance, if all ki = 0, then Vj-i = and the 
equation ()4.19p reduces to Vj^i = for all j = 1, . . . , ttiq. Hence we obtain that in this 
case the condition is that the coefficients in front of {z — w^)~^ and {z — w^) in the 
expansion of v(z) in Laurent series in {z — w^) should vanish for all j = 1, . . . , uiq. 

Thus, a connection dz+u{z) given by formula (j6.16p with the numbers w^, wj satisfy- 
ing the Bethe Ansatz equations (|6.2Up - (|6.2ip gives rise to an operator (|6.17p such that 
its coefficients in the expansion at w^,j = 1, . . . ,mo, satisfy the equations (|6.19p and 
(j4.19p . (However, we do not expect the converse to be true in general.) Our proposal is 
that the generic common eigenvalues of the affine Gaudin Hamiltonians on the tensor 
product (S>iLi of Verma modules over 5(2 are encoded by 5[2-opers of this kind. In 
general, some degenerations of these opers will also correspond to some eigenvectors, 
as in the shift of argument Gaudin model (see Section [4.4p . 

It is useful to express the operator (|6.17p as a second order differential operator 
(compare with Section Applying the (formal) gauge transformation by the element 
nili(-z - Zi)^^"^, we obtain 

N 



which may be rewritten as the following second order operator with spectral parameter 

N 

dl - V{Z) - \\{{Z - Zif^ . 
i=l 

Alternatively, we may rewrite (j6.17p as the following partial differential operator: 

{dz + Y.^^,^d^ -v{z)-\. 
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6.10. Integrable modules. Next, we consider the common eigenvalues of the affine 
Gaudin Hamiltonians on the tensor product of integrable g-modules. Recall that for 
each integral dominant weight v we have an integrable module V^, which is the irre- 
ducible quotient of the Verma module My. By analogy with the finite-dimensional case, 
we propose the following description of the common eigenvalues of the affine Gaudin 
Hamiltonians on the tensor product '^Vi- 

Conjecture 6. There is a bijection between the set of common eigenvalues of the 
affine Gaudin Hamiltonians on the tensor product (^^i V^- of integrable modules over 
Q (counted without multiplicities) and the ^g-opers on which are Vi-regular at Zi,i = 
1, . . . ,N, and have regular singularity at oo. 

Note that such opers necessarily have trivial monodromy around each of the points 
Zi,i = 1, . . . ,N, and hence they have no monodromy on P"*^. 

If the Bethe Ansatz equations ()6.4p are satisfied and the projection of the Bethe 
vector (16. 3p onto V^. is non-zero, then we obtain an eigenvector of the affine 

Gaudin Hamiltonians in (^^^ V^.. A necessary condition for this is that its weight 

N m 

(6.22) 7 = E^^-E"^. 

1=1 3=1 

is dominant. The corresponding eigenvalue should then be encoded by the ^g-oper 
obtained by applying the Miura transformation to the Cartan connection (|6.15p . As 
shown in [F2], such an oper satisfies the conditions of Conjecture [6j This provides the 
first consistency check for this conjecture. 

Note that we allow the possibility that there exist eigenvectors other than the Bethe 
vectors in (^^^ Vy- . In this case we still believe that the common eigenvalues of the 
affine Gaudin Hamiltonians are encoded by the ^g-opers as in Conjecture [H However, 
the eigenvalues corresponding to the non-Bethe eigenvectors should be encoded by 
opers that cannot be expressed as the Miura transformation of a connection of the 
form (j6.15p , of the kind described at the end of Section 14.41 

Finally, if the weight (I6.22P is not dominant, then the corresponding Bethe vector 
in Vy. is zero. However, it is still interesting to consider the corresponding space 

of solutions of the Bethe Ansatz equations (16. 4[) . As shown in [F2], Theorem 5.7, for 
each fixed weight 7 this set may be identified with a disjoint union of open subsets of 
a finite-dimensional ^i?_-orbit (Schubert cell) in the affine fiag variety ^G/^B^. We 
will consider a special case of this correspondence in the next subsection. 

6.11. Special case: rational solutions of the KdV hierarchy. Consider the spe- 
cial case of the above construction when g = 5(2, = 1, the highest weight vi = 0, 
and the level ki = 0. According to formula (j6.18p . the corresponding generic sl2-opers 
may be represented by the following second order differential operators: 

mo 2 

(6-23) ^'-Et 0^-^' 

^ iz - w^'Y 

j=l 3 
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which satisfy the condition 

(6.24) ^^^^ = 0, j = l,...,mo. 

These equations (analogous to the Bethe Ansatz equations (|6.4|) ) mean that the (z—w^)- 
coefficient Vj,i of the expansion oiv{z) at z = w'j vanishes, which is equivalent to (j4.19p 
in this case. This condition means that the operator (j6.23p has no monodromy on 
for all values of the spectral parameter A. 

The operators (j6.23p were introduced by Airault-McKean-Moser in |AMMj as the 
generic rational solutions of the KdV hierarchy. They were also encountered by Duis- 
termaat and Griinbaum |DG| in their study of the bispectral problem for second order 
differential operators. Here we have obtained these operators in a different way; namely, 
in the framework of the affine Gaudin model in the special case when N = = 
and ki = 0. Following the argument of |F2j . we obtain that the set of these operators 
may be identified with an open subset of the 5'L2[[A]]-orbit 



in the affine Grassmannian S L2{{X)) / S L2[[X]]. 

The points outside this locus correspond to the sl2-opers in which some of the points 
coalesce. They are represented by the more general second order operators of the 
form 

(6-25) d! - f: f^^L+il _ A, 

j=l 3' 

where € Z>o, which satisfy the no-monodromy condition 

(6.26) E iJo^^'m,+i = j = l,...,mo;p, = l,...,£,-. 

These are precisely the operators constructed in |DG] (see formulas (3.31) and (3.32)), 
where it was shown that they correspond to the most general rational solutions of the 
KdV hierarchy decaying at oo. The fact that the set of all operators of this form is 
in bijection with the affine Grassmannian S'L2((A))/S'L2[[A]] is in agreement with the 
well-known description of the polynomial tau-functions of the KdV hierarchy as points 
of the affine Grassmannian 

The rational solutions of the g-KdV hierarchies for other affine Kac-Moody algebras 
may be interpreted in a similar way. 

7. Concluding remarks 

Let us summarize what we have learned so far. Our point of departure was the 
observation that soliton systems, such as the KdV and AKNS hierarchies, are special 
cases of the Gaudin models associated to affine Kac-Moody algebras. Therefore the 



"'^'''An identification of tlie space of rational solutions of the KdV hierarchy with the space of afhne 



opers of this kind may also be obtained from the results of [BF| . 



64 



BORIS FEIGIN AND EDWARD FRENKEL 



quantization of these models should be understood in this framework. In order to better 
understand these models, we look at their simplified versions; namely, the Gaudin 
models associated to finite-dimensional simple Lie algebras. We then use this analogy 
to describe various elements of the affine Gaudin models, and in particular the quantum 
soliton systems. In this paper we have made the first steps in this direction. 

Let us first discuss the classical soliton systems. These are integrable systems defined 
on the space of connections 

dt + A{t), A{t)eg{{t)). 
We apply a shift of argument and consider the space of connections of the form 
(7.1) dt + A{t) + zx, A{t)eg{{t)). 

The shift parameter z plays the role of a new variable, and so these integrable systems 
are properly understood as integrable systems on the double loop algebra in variables 
t and z, as was already noted in jRSlj . The classical commuting Hamiltonians are 
obtained from the functions (p{M{z)), where M{z) is the monodromy matrix of (j7.ip 
and 93 is a function on G invariant under conjugation. Expanding these functions in 
z, we obtain local and non-local Hamiltonians. It is important to realize that, when 
viewed as a function on the entire double loop algebra, ip{M{z)) is actually central, 
with respect to the natural Kirillov-Kostant structure. 

How can we quantize this system? It is clear that the quantum Hamiltonians should 
be obtained from quantum versions of the functions (/?(M(z)), which should be in the 
center of a proper completion of the enveloping algebra of the double loop algebra. 
Therefore to describe the spectrum of these Hamiltonians we need to understand the 
structure of this center. 

At this point it is instructive to look more closely at the "baby version" of this 
model; namely, the Gaudin model associated to a finite-dimensional simple Lie algebra 
0. In this case the naive phase space is g*, but because we make the shift of argument, 
the true phase space is the dual of the loop algebra q{{z)). The classical Hamiltonians 
appear from the z-expansion of the central functions on g((z))*, obtained from invariant 
polynomials on q* . Naively, we expect that the quantum versions of these functions are 
central elements in a completed enveloping algebra of q([z)). But there is an important 
twist to this story: we need to replace q{{z)) by its central extension; its completed 
enveloping algebra then has a large center at the critical level. 

This center turns out to be isomorphic to the algebra of functions on the space 
OpLg{D^ ) of ^g-opers on the punctured disc, where is the Langlands dual Lie algebra 
of g |FF3l [F6j . This is a non-trivial result, which is closely related to the geometric 
Langlands correspondence |BDll [F5] . This implies [F31 IFFT] that the spectrum of the 
algebra of quantum Hamiltonians of the Gaudin models associated to g is isomorphic 
to the space of ^g-opers on with prescribed singularities at finitely many points 
(depending on the particular model). 

Now, we conjecture that the algebras of quantum Hamiltonians of the affine Gaudin 
models exhibit similar features. Namely, their spectra should be expressed in terms of 
^g-opers on with prescribed singularities at finitely many points, where ^g is the 
Langlands dual Lie algebra of g. These ^g-opers may be described more concretely as 
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differential operators. For example, for q = 5(2 these are Schrodinger operators with 
spectral parameter. 

We remark that in this picture opers appear in two different contexts, and this 
may cause some confusion. First of all, the classical fl-KdV system has the space of 
g-opers (or g-opers) on the punctured disc as the phase space. Then we quantize it 
and look at the spectra of the corresponding algebra of quantum KdV Hamiltonians. 
We conjecture that those are encoded by ^g-opers on P^. So opers appear again, but 
they have nothing to do with the g-opers entering the definition of the classical KdV 
system! For one thing, they correspond to an a priori different Lie algebra; namely, the 
Langlands dual Lie algebra ^g. Besides, they are defined on (on which the spectral 
parameter z is a global coordinate), not on a punctured disc. 

In this paper we have presented explicit conjectures for the ^g-opers encoding the 
spectra of the quantum Hamiltonians of the affine Gaudin models associated to g. We 
have checked that our proposal coincides with the proposal of |BLZ5j in the case of the 
quantum KdV Hamiltonians. 

What can we learn from these conjectures? First of all, placing the quantum soliton 
systems in the context of affine Gaudin models is useful, because we may apply the tools 
of the Gaudin models to investigate the soliton systems. An example of such a tool is 
Bethe Ansatz, which we use here to produce explicit formulas for the eigenvectors of 
the quantum Hamiltonians. Second, this points us to a conjectural description of the 
spectra of the quantum Hamiltonians in terms of opers, which is consistent with and 
generalizes the proposals obtained by other methods. Finally, we place the quantum 
soliton systems in the context of Langlands duality. We hope that this will allow us to 
understand better the true meaning of quantum soliton system and their relations to 
other areas of mathematics and physics. 



8. Appendix. Proof of Lemma [2] 



For any simple Lie algebra g we have the Segal-Sugawara current 



S{w) = ^ ^ {■.ea{w)fa{w): + :fa{w)ea{w):) + ^^:hi{w)h\w): 




where {/i*} is the dual basis to the coroot basis {hi} of f). Let us set 




aGA+ 



and 




i=i 



Now denote by 



Sa{w) = ^ Sa,n'W 



-n-2 
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the element S{w) associated to the 3(2 subalgebra fla C 5 generated by and fa, and 
by 

the element S{w) associated to the rank two Lie subalgebra Qai3 C g whose root system 
Aq,^ C a is spanned by the roots a,/3 € A_(_. We will use a similar notation for the 
corresponding currents S{w) and S{w). 

Adapting the argument proving the commutativity of the DMT Hamiltonians in the 
finite-dimensional case (see [TL] ) , we reduce the statement of part (1) of Lemma [2] to 
checking the following identity: 

(8.1) [S-y,o,Sap,o]=0, V7GA+^. 
We have 

5'7,o = 'S'7,0 — 5*7,0 — 70) 

Sal3,0 = Sapfl — Saf3,0 — 2Pq/3,0- 

Here 7 (resp., p) denotes the element of f) corresponding to the root 7 (resp., the half- 
sum of positive roots of Qais) in f)* under the identification f) ~ {)* corresponding to hq. 
We write 70 (resp., po) for 7 1 (resp., p iX) 1) in q. 

Since all summands in the right hand sides of these formulas commute with f) 1 C g, 
and so in particular with Pa/3,0, we obtain that (18. ip is equivalent to 

[5*7,0) SafS^O — 5q,/3,o] = 0. 

Actually, we will show that 

[Syfi, Saf3,o\ = [5*7,0, 5'q,/3,o] = 0. 

By using the formalism of operator product expansion, it is easy to derive these iden- 
tities from the following identities: 

(8.2) Sai3,~iS^,-2Vk G ImT, 

(8.3) Saf3-iS^-2Vk G ImT 

in the vertex algebra Vfc(fl), where is the vacuum vector and T is the translation 
operator, in the notation of [FBj . (The sought-after commutation relations will follow 
after we make the change of variables if = e" — 1 in the corresponding operator product 
expansion.) 

To check ()8.ip . recall that the Segal-Sugawara operator Sap^-i acts on the completed 
enveloping algebra of g^^ as a multiple of the derivation T = —tdt- This implies formula 

(j8.2[) . To prove formula (18. 3p . note that S^-2Vk = e^-if^-iVk and Sap-i is a linear 
combination of 



(8.4) 



7-i-n7n and 
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where /i G f) is orthogonal to 7. Therefore the second term, apphed to S^^-2Vk, gives 
0, and the first term gives 

^7_i_n7n • ey-if^-iVk = 27_i7_2Ufc = T ■ {"y^ifvk. 

This imphes formula (j8.3p and hence proves part (1) of Lemma [2j 

The commutativity of the operators given by formula (I3.12p follows from the 
corresponding analogue of the identity (|8.1|) . which is proved in the same way as above. 
In order to complete the proof of part (2) of Lemma O it remains to show that the 
operators commute with \). This, in turn, follows from the formula 

which is verified in a straightforward fashion. 
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